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ABSTRACT 

We present UaaeUeoeA 25 L^i 60 follow-up HST observations of 16 z ^ 3 candidate LyC emitters in 
the HS1549-I-I919 field. With these data, we obtain high spatial-resolution photometric redshifts of all 
sub-arcsecond components of the LyC candidates in order to eliminate foreground contamination and 
identify robust candidates for leaking LyC emission. Of the 16 candidates, we find one object with a 
robust LyC detection that is not due to foreground contamination. This object (MD5) resolves into 
two components; we refer to the LyC-emitting component as MD5b. MD5b has an observed 1500A 
to 900A flux-density ratio of {Fuv/FLyc)obs = 4.0 ± 2.0, compatible with predictions from stellar 
population synthesis models. Assuming minimal IGM absorption, this ratio corresponds to a relative 
(absolute) escape fraction of = 75 — 100% {f^fabs = 14 — 19%). The stellar population 

fit to MD5b indicates an age of < 50Myr, which is in the youngest 10% of the FIST sample and 
the youngest third of typical z ^ 3 Lyman break galaxies, and may be a contributing factor to its 
LyC detection. We obtain a revised, contamination-free estimate for the comoving specific ionizing 
emissivity at z = 2.85, indicating (with large uncertainties) that star-forming galaxies provide roughly 
the same contribution as QSOs to the ionizing background at this redshift. Our results show that 
foreground contamination prevents ground-based LyC studies from obtaining a full understanding of 
LyC emission from z ~ 3 star-forming galaxies. Future progress in direct LyC searches is contingent 
upon the elimination of foreground contaminants through high spatial-resolution observations, and 
upon acquisition of sufficiently deep LyC imaging to probe ionizing radiation in high-redshift galaxies. 
Subject headings: galaxies: high-redshift - intergalactic medium - cosmology: observations - diffuse 
radiation 


1. INTRODUCTION 

The sources responsible for cosmic reionization are 
still not well understood. Evidence that quasars 
QSOs) cannot be solely responsib l e for reionization 
Fontanot et al.l . 120121 : iClikman et al.l . 120111 : iSiana et al.l . 
l2008f l has prompted many searches for ionizing Ly¬ 
man continuum (LyC) emission from star-forming galax¬ 
ies. While the IGM at z > 6 is opaque to 

LyC photons and prevents direct observations of LyC- 
emitting galaxies during the epoch of reionization, 
many studies have attempted to detect lower-redshift 
analogs to galaxies responsible for reionization. Al¬ 
though IGM transmission is high est in the loc a l uni¬ 


verse, studies at z < 2 (e.g., IGrimes et al.l. 12007 

200S 

:ICowie et al.l. 120091: iBridee et al.l.l2010l:ISiana et al. 

2007 

. |2010D have yielded yery few detections of LyC 


emission, wi t h on l y three objects identif i ed to date 
(jLeitet et al.l . 1201 iL I20I3I : iBorthakur et al.l . 1201411 . At 
redshift z 3 — 4, the search for LyC-emitting galax¬ 
ies has appeared to be more fruitful. However, even 
though the e xamination o f hundreds of galax i es (in 
works such as iSteidel et al.L 12001: iShaolev et 12006 : 
llwata et'ffil . 120091: iNestor et al.l . (20131: iMostardi et ah . 
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2f)Ta ISiana et aLL iVanzella, et all. 1201 OH IMl 

2015f ) has yielded many promising LyC-emitting can- 
didates, there exist on ly two robust detections 
(IVanzella et al.l. 

Amassing large samples of LyC detections in high- 
redshift star-forming galaxies has been difficult for sev¬ 
eral reasons. First, large parent samples of high-redshift 
galaxies must be identified and confirmed spectroscopi¬ 
cally, requiring extensive galaxy surveys (often ground- 
based) and time-consuming spectroscopic follow-up. Sec¬ 
ond, it is necessary to probe flux in the LyC spectral 
region for these galaxies, either with deep spectroscopy, 
also very time-intensive, or through narrowband imag¬ 
ing in a filter just blueward of the Lyman limit, in which 
it is difficult to match a single narrowband filter to the 
LyC region for many galaxies at once. Even after poten¬ 
tial high-redshift LyC-emitting candidates are identified, 
there remains the possibility that apparent LyC emission 
is actually due to a lower-redshift interloper along the line 
of sight, which cannot be distinguished in ground-based, 
seeing-limited data. 

One method that has proven successful at identify¬ 
ing potential LyC-emitting galaxies is narrowband LyC 
imaging of galaxy protoclusters. Large ground-based sur- 
veys of UV-selected star-forming galaxies at z ^ 2 — 3 
(|Steidel et al.l . 1200^ 120041201 lUR.eddv et al.l . 120081 1 have 
identified and spectroscopically confirmed thousands of 
high-redshift star-forming galaxies. These surveys have 
also located galaxy protoclusters, areas on the sky with 
large overdensities of galaxies at similar redshift. A 
very effective way to simultaneously probe the LyC of 
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large samples of galaxies at the same redshift is to per¬ 
form deep imaging through a narrowband filter tuned 
to the Ly C spectral r egioii at the protocluster red- 
shift (e.g.. Ilwa,ta et al.l . 120091 : iNestor et al.l . I2011L 1201,11 
iMostardi et alJ . 1201, 11 1. 

Initially, these protocluster studies were entirely based 
upon ground-based data with seeing FWHMs of 0'.'7 — 
I'/O, and thus suffered from the possibility of fore¬ 
ground contamination. iVanzella et al.l (l2010all demon¬ 
strated that statistical simulations modeling the distri¬ 
bution of foreground galaxies result in high rates of 
foreground contamination for high-redshift objects in 
ground-based studies. While simulations can account for 
co ntamination s t atisti c ally in LyC-emi t ting s amples (as 
in INestor et al.l . 120111: IMostardi et al.L 120131) . contami¬ 
nants cannot be eradicated on an individual basis. As 
two of the main goals of LyC studies are to determine 
the mechanism of LyC photon escape from the interstel¬ 
lar medium (ISM), and to identify additional features of 
LyC-emitting galaxies that may enable their identifica¬ 
tion through other means, it is crucial to identify robust, 
individual candidates for LyC emission where foreground 
contamination has been ruled out. 

Eliminating contaminants is a complex process. High- 
resolution imaging shows that the majority of high- 
redshift galaxies are not morphologically simple, but 
are composed of m ultiple compact c lumps and/or dif¬ 
fuse emission (e.g., iLaw et all . I2007D . Contamination 
can only be firmly ruled out if the redshifts of indi¬ 
vidual galaxy clumps are measured, and if the clump 
associated with LyC emission is confirmed to be at 
the redshift of the target galaxy. In order to address 
the issue of contamination in the narrowb and LyC sur¬ 
vey of the z = 3.1 SSA 2 2a pr otocluster (jNestor et al.l . 
l2?)ril IWl . iSiana etldl (|201,5D obtained near-IR spec¬ 
troscopy with Keck/NIRSPEC to measure the spectro¬ 
scopic redshifts of the sub-arcsecond components of 5 
LyC candidates. These authors found two foreground 
contaminants, one galaxy with a misidentified redshift, 
and two galaxies that could not be definitively con- 
hrmed as LyC-emitte rs. For galaxies at sli ghtly higher 
redshifts (z ~ 3.7), IVanzella et al.l (|2ni2[ ) used pho¬ 
tometric redshifts obtained throug h the high-resolution 
imaging in the CANDE LS survey (iGrogin et al.l . 120111 : 
iKoekemoer et al.l . 120111) to analyze the sub-arcsecond 
clumps of 19 candidate LyC-emitters, finding 18 contam¬ 
inants and one bona-fide LyC emitter. These two studies 
have shown that both methods — high-resolution spec¬ 
troscopy and photometric redshifts — are effective ways 
to lo cate foreground c ontam inants. 

In IMostardi et al.l (I2013D (hereafter M13), we per¬ 
formed a narrowband LyC imaging surve y of a galaxy 
proto cluster in the HS1549-I-1919 field (jSteidel et al.l . 
I 2 OIII) . In this work, we present follow-up observations 
with HST for the sample of candidate LyC emitters 
in M13 with the goal of using photometric redshifts to 
eliminate contaminants from the LyC emitter sample. 
The HS1549 protocluster has a redshift-space overden¬ 
sity of Sgai 5 at z = 2.85 ± 0.03, and this “spike” 
redshift coincides with t hat of a hyperluminous QSO 
(|Trainor fc Steidel I2012D . More than 350 UV-selected 
galaxies have been identified in the HS1549 field, ~ 160 
of which have been spectroscopically confirmed at 1.5 < 


z < 3.5. Additionally, narrowband imaging with a 4670A 
filter tuned to the wavelength of Lya at the redshift spike 
has revealed ^300 p otential Lya Emitters (L AEs) and 
several Lya “blobs” (|Steidel et al.l . l200'ol 120111) . In M13, 
we used a narrowband filter (NB3420) tuned to wave¬ 
lengths just below the Lyman limit at z > 2.82, thus 
observing the LyC spectral region for hundreds of Ly¬ 
man break galaxies (LBGs) and LAEs at z > 2.82, in¬ 
cluding 49 LBGs and 91 LAEs with spectroscopic con¬ 
firmation. We identified putative LyG emitters in the 
NB3420 imaging, and also performed a stacking analysis 
of objects undetected in the NB3420 filter (measuring no 
signal). Although we found an NB3420 detection rate of 
~ 9% in both the LBG and LAE samples, simulations 
indicated that 40-75% of the individual NB3420 detec¬ 
tions may have resulted from foreground contamination, 
highlighting the need for further work to disentangle true 
LyG emitters from low-redshift contaminants. 

Our aims in the current work are two-fold. Our pri¬ 
mary goal is to address the question of foreground con- 
ta mination in the M 13 LyC-emitter sample. Ideally, as 
in ISiana et "mI (|2015l) . we would obtain spatially-resolved 
spectroscopy in the vicinity of each putative LyC detec¬ 
tion, with resolution of < 0'.'5, in order to measure the 
redshifts of all components. However, ground-based op¬ 
tical spectroscopy probing the rest-frame UV provides 
insufficient spatial resolution, and ground-based near- 
IR spectroscopy of rest-frame optical nebular emission 
lines (with or without the assistance of adaptive optics) 
is not feasible because at z ~ 2.85 the strongest nebular 
emission lines (Ha, [OHIJA5007, H/3, and [OHJA3727) are 
lost either in the thermal background or gaps in atmo¬ 
spheric transmission. Therefore, we have identified high- 
resolution, multi-band HST imaging as the best method 
for estimating spatially resolved photometric redshifts for 
the individual galaxy subcomponents associated with ap¬ 
parent LyC emission. Our second goal is to analyze the 
properties of galaxies we have verified to be true sources 
of LyC emission (such as their morphologies, stellar pop¬ 
ulations, and the ratio of their ionizing to non-ionizing 
flux densities) with respect to properties of star-forming 
galaxies without LyC detections. Such an analysis will 
help determine whether galaxies with high escape frac¬ 
tions of ionizing photons have different intrinsic prop¬ 
erties from those of galaxies without detectable leaking 
ionizing radiation, and may provide insight into star for¬ 
mation and the structure of the ISM in high-redshift 
galaxies. 

In addition to presenting high-resolution, multiwave¬ 
length follow-up HST observations of high-redshift LyC- 
emitting candidates at z ~ 2.85 from M13, we discuss the 
implications for continuing searches for ionizing radiation 
in star-forming galaxies. The paper is organized as fol¬ 
lows. Section [2] describes our methodology, while the 
galaxy sample and HST observations are presented in 
Section[3] In Section|3]we describe the techniques used to 
reduce the HST imaging data and perform photometry 
on the sub-arcsecond components of each galaxy. Section 
[5] contains a discussion of our methods of fitting photo¬ 
metric redshifts to the data, as well as an analysis of each 
candidate LyC emitter and any sources of foreground 
contamination. In Section m we analyze the properties 
of our best candidate for true LyC emission (MD5) with 
respect to the larger population of star-forming galax- 
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Figure 1. C/ 336 , Veoe, A 25 , and iiCieo filters superimposed 

over a model 2 = 2.85 galaxy spectrum. The model spectrum 
is a solar- metallicity 100 Myr const ant star-formation rate 
model from IBruzual fc Charlotl ||200^ . reddened to E{B — V) 
= 0.15 with the lCalzetti et aH (120001 1 attenuation curve. The 
C /336 filter probes the LyC spectral region, but with some 
contamination (~ 20%) redwards of the Lyman limit. The 
Veoe filter probes the rest-frame non-ionizing UV continuum, 
and the J 125 and iiCieo filters probe flux on either side of the 
Balmer break. With both the Lyman and Balmer breaks 
sensitively probed at 2 = 2.85, we can determine photometric 
redshifts for individual sub-arcsecond components of galax¬ 
ies within our LAE and LBG samples. These photometric 
redshifts enable us to determine whether the NB3420 detec¬ 
tions associated with these galaxies are true LyC emission, or 
contamination from foreground galaxies, 
ies. In Section [71 we consider the broader implications of 
these results, and the prospects for future direct searches 
for LyC radiation at high redshift. We summarize our 
results in Section |8l Throughout the paper we employ 
the AB magnitude system and assume a cosmology with 
= 0.3, Ha = 0.7, and Hq = 70 km s“^ Mpc“^. At 
z = 2.85, 1"corresponds to 7.8 proper kpc. 

2. METHODOLOGY 

To evaluate the amount of foreground contamina¬ 
tion in M13 using photometric redshifts, we selected 
four HST filters (WFC3/UVIS C/ 336 , ACS/WFC Feoe, 
and WFC3/IR J 125 and LCieo) designed to probe the 
strengths of the Lyman and Balmer breaks at z 
2.85. Figure [B shows the locations of the C/ 336 , 
Vsns, A 3 . 6 , and H^m filte rs superimposed over a model 
iBruziial fc Charlc^ (j2003[l galaxy spectrum redshifted to 
z = 2.85. The choice of J 125 and .H 160 filters is particu¬ 
larly powerful for this test, due to the observed wave¬ 
length of the Balmer break at z = 2.85 {Xobs,BB = 
1.4/im). This wavelength corresponds exactly to the red 
cut-off of J 125 and the blue cut-off of i7i6o- Therefore, the 
Ji 25 — Hiqo color is very sensitive to the presence of the 
Balmer break, and provides information about the age 
of the stellar population. Additionally, the Veos — ^125 
color probes the rest-frame UV slope at z ^ 2.85, provid¬ 
ing information about the stellar populations and dust 
extinctio n. At z < 2, where m ost of the contaminants 
found by iVanzella et al.l (|2012h are located, the J 125 fil¬ 
ter falls entirely on the red side of the Balmer break, 
and therefore, extremely flat J 125 — 7Ci60 colors are ex¬ 
pected, with the Balmer break falling instead between 
Vqog and J 125 . At the other end of the spectrum, the 
C /336 filter does not lie entirely bluewards of the Lyman 


limit at z = 2.85 (only 80% of its wavelength range falls 
below the Lyman limit) and thus does not exclusively 
probe the LyC spectral region. However, given that 
the Lyman break passes through this filter at redshifts 
z = 2.40 — 2.95, the C /336 — ^606 color should be sensitive 
to the magnitude of the Lyman break at z = 2.85, and 
to its absence at significantly lower redshift. With both 
Lyman and Balmer breaks sensitively probed at z = 2.85 
using C /336 V606</i257/i 60 photometry, we can distinguish 
between true sources of LyC emission, and those NB3420 
detections attributable to non-ionizing radiation at lower 
redshift. 

Several complications may arise from the method of es¬ 
timating the photometric redshifts using only these four 
filters. First, galaxies with leaking LyC emission may 
have intrinsically high ratios of escaping ionizing to non¬ 
ionizing radiation, and may not exhibit as strong a Ly¬ 
man break as expected from normal z = 2.85 galaxies. 
We keep this caveat in mind during our subsequent anal¬ 
ysis, with the understanding that the C /336 ~ ^606 color 
may not be accurately represented by the models. Sec¬ 
ond, young galaxies (< 20 Myr) may not have a signif¬ 
icant Balmer break. The combination of these two sce¬ 
narios may result in a young, high-redshift, LyC emitter 
with a relatively featureless spectrum, making it difficult 
to distinguish such galaxies from low-redshift contami¬ 
nants. In our photometric redshift analysis, we highlight 
cases of galaxies with ambiguous SED shapes and the 
possible stellar populations these SEDs may indicate. 

3. SAMPLE AND OBSERVATIONS 

The HS1549 galaxy sample discussed in M13 consists of 
49 LBGs and 91 LAEs, all spectroscopically confirmed at 
z > 2.82. There are 9 objects part of both the LAE and 
LBG samples; for simplicity in displaying our data, we 
group these objects with the LBGs0 The redshift limit 
of z > 2.82 ensures that the NB3420 filter is sensitive 
only to LyC emission, with no contamination from flux 
redward of the Lyman limit. Out of these galaxies, 5 
LBGs and 7 LAEs have NB3420 detections within 1'.'9 of 
the non-ionizing UV (~ 1500A) centroid of the galaxy, 
indicating potential LyG emission if there is no lower- 
redshift galaxy contaminant along the line of sight. 

In M13 we also present an additional sample of 33 
photometric LAE candidates (no spectroscopic confirma¬ 
tion) in the HS1549 field whose magnitudes in the nar¬ 
rowband filter used to select LAEs (NB4670) are in the 
same range as those of the spectroscopically-confirmed 
LAEfl While these objects are not included in the anal¬ 
ysis of M13 because the lack of spectroscopic redshifts in¬ 
creases the possibility of contamination, 10 of them are 
reported to have NB3420 detections. 

Einally, in M13 we present 8 spectroscopically con¬ 
firmed LAEs with NB3420 detections that are either 
fainter than the magnitude limit of the main LAE sam¬ 
ple (to 467 o > 26; 5 objects) or were selected by their 
G—NB4670 color rather than U—NB4670 (3 objects). 

® One of the objects in both the LAE and LBG samples is the 
putative LyC-emitter MD12/lae3540y which is discussed in detail 
in M13. In this work, we refer to this object simply as MD12. 

^ All LAEs within the spectroscopic LAE sample of M13 had 
^4670 ^ 26. Thus, the photometric LAE sample of M13 was de¬ 
fined to be any additional LAEs with m 467 o < 26 that did not have 
spectroscopic confirmation. 
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Figure 2. Footprints of HST pointings superimposed on the 5' x 1' LRIS NB3420 (LyC) image. Blue squares indicate the 
footprint of t/ 336 , green squares indicate the footprint of Veoe, and red squares indicate the footprints of J 125 and Cyan 

and yellow circles indicate the positions of LBGs and LAEs with NB3420 detections that lie within the footprints of the HST 
images. The naming scheme for LBGs is presented in ISteidel et al.l (1200311 . 
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Table 1 

Description of Samples 


Sample® 

Ntot'^ 

^NB3420 WuvjhT 

^NB3420 ^^uvjh]’^ 

LBGs (with Zspec) 

49 

5 [4] 

44 [12] 

LAEs (with Zspec) 

82® 

7 [6] 

75 [18] 

LAEs (without Zspec) 

33 

10 [4] 

23 [7] 

Faint LAEs (with Zspec) 

— 

8 [2] 

-- 

GNBs (with Zspec) 

— 

3 [0] 

-■ 


The samples described in M13. 

^ The total number of galaxies in the M13 sample. 

^ The number of galaxies with NB3420 detections. In brackets, we indicate the 
number of these galaxies for which we have obtained t/336V' 606 ' 7 i 25 ^f 160 imaging. 

The number of galaxies without NB3420 detections. In brackets, we indicate the 
number of these galaxies for which we have obtained t/336V606'/i25//i60 imaging. 
® The number of LAEs that are not part of the LBG sample, i.e., 91 LAEs minus 
the 9 overlap objects (which are listed here as part of the LBG sample). 


We did not include these objects in the analysis of Ml3 
because we had not assembled a complete and unbiased 
spectroscopic sample of LAEs with TO 4670 > 26 or LAEs 
selected by their G—NB4670 color. 

In the follow-up observations presented in this work, 
we considered all objects with NB3420 detections pre¬ 
sented in M13 and chose HST pointings oriented to 
maximize the number of these potential LyC-emitting 
targets on the image footprints. As the fields of view 
attained by the various HST instruments employed for 


these observations (ACS/WFC: 315x315; WFC3/UVIS: 
219x217; WFC3/IR: 213x211) are much smaller than the 
Keck/LRIS field of view (5'x7'), we were unable to ac¬ 
quire imaging for all of the potential LyC-emitters. How¬ 
ever, with two separate HST pointings in each of the 
4 filters (see Figure [2]), we obtained 1 / 336 V 606 '^i 2577 i 60 
photometry for all but two of the spectroscopically con¬ 
firmed galaxies in the main sample: 4/5 LBGs and 6/7 
LAEs. We also obtained 1 / 336 Veoe/i 257 /i 60 photometry 
for 4/10 photometric LAE candidates, and 2/6 of spec- 
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Table 2 

HST Imaging Observations 


Filter 

A)' 

PSF FWHM 
(") 

Depth® 

(mag) 

Exposure 

(s) 

Pixel Scale*’ 
("/pixel) 

Pixfrac*’ 

U 336 

3355 

0.081 

29.20 

14176 

0.025 

0.7 

Veoe 

5921 

0.092 

30.22 

11848 

0.03 

0.7 

/125 

12486 

0.178 

28.93 

7835 

0.075 

0.7 

Hwo 

15369 

0.186 

28.63 

7835 

0.075 

0.7 


^ The Sct limiting depth obtained in a circular aperture with a diameter of 1.5 
times the PSF FWHM. 

^ The AstroDrizzle parameters for pixel scale (plxscale) and pixf rac listed here 
represent the parameters used to attain optimal resolution. 


NB3420 V F336W F606W F125W F160W UVJ Pixel Map 



Figure 3. 5” X 5" postage stamp images of spectroscopically-confirmed LBGs with NB3420 detections and imaging in all 

four HST filters. From left to right, objects are displayed in the LRIS NB4670—F continuum-subtracted image (indicating 
Lya emission; MD12 only), LRIS NB3420 (LyC emission), LRIS V (non-ionizing UV continuum), HST U 33 Q (a combination 
of LyC and non-ionizing UV), HST Veoe (non-ionizing UV continuum), HST J 125 (optical, bluewards of the Balmer break), 
and HST iLieo (optical, redwards of the Balmer break). The penultimate column shows a color-composite image of HST 
U 336 (blue), Veoe (green), and J 125 (red). The final column shows SExtractor segmentation maps of the pixels used for the 
photometry of each labeled galaxy component; the arbitrary color scale indicates component edges when two components are 
adjacent. Only pixel maps for the labeled objects are shown. Red circles (I'.'O diameter) indicate the centroid of the NB3420 
emission, and blue circles indicate the centroid of the Lya emission for MD12. Photometry was performed individually on 
sub-arcsecond components associated with each LBG, and all components are labeled in the Veoe image. Postage stamps follow 
the conventional orientation, with north up and east to the left. 


troscopically confirmed LAEs with m 467 o > 26, total¬ 
ing 16 galaxies with NB3420 detections covered in all 
four filters. Eight additional objects with NB3420 de¬ 
tections are covered by at least one HST filter (usually 
Vqog, which has the largest field of view), and may be 
examined morphologically, although it is not possible to 
obtain photometric redshifts for galaxies without imag¬ 
ing in all four filters. These objects are presented in 
Appendix B. In addition to the 16 galaxies with NB3420 
detections, we acquired C/ 336 V 606 <^i 25 .ffi 60 imaging for 30 
spectroscopically confirmed galaxies at z > 2.82 without 
NB3420 detections (12 LBGs and 18 LAEs), which allows 
us to calibrate our photometric redshift fitting methods 
on galaxies without LyC detections and facilitates the 
differential analysis of the stellar populations of galaxies 
with and without LyC detections. Finally, 50 additional 


spectroscopically confirmed galaxies at z > 2.82 without 
NB3420 detections (30 LBGs and 20 LAEs) were par¬ 
tially covered by our suite of HST imaging. Table [T] 
summarizes the HST coverage of the samples. 

In total, we obtained 5 orbits for each of the 
WFC3/UVIS C /336 and ACS/WFC Veoe pointings and 
3 orbits for each of the WFC3/IR J 125 and //leo point¬ 
ings as part of HST Program ID 12959 (PI: A. Shap- 
ley) between 2012 December and 2013 August. Table 
[2] lists details of the observations. Individual exposures 
were half-orbit ('^IdOO seconds) for C /336 and quarter- 
orbit (^600 seconds) for Veosj >/i 25 , and LCieoj and total 
exposure times per pointing were 14 ks (C/ 336 ), 12 ks 
(heoe)) and 8 ks (J 125 & /Ci 6 o)- We used a combination 
of the WFC3/UVIS DITHER-BOX and DITHER-LINE 
patterns for C/ 336 , the ACS/WFC DITHER-BOX pat- 
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Figure 4. 5” X 5” postage stamp images of spectroscopically-confirmed LAEs with NB3420 detections and imaging in all 

four HST filters. From left to right, objects are displayed in LRIS NB3420 (LyC emission), LRIS NB4670—E (indicating 
Lya emission), LRIS V (non-ionizing UV continuum), HST C /336 (a combination of LyC and non-ionizing UV), HST Veoe 
(non-ionizing UV continuum), HST J 125 (optical, bluewards of the Balmer break), and HST LCieo (optical, redwards of the 
Balmer break). The penultimate column shows a color-composite image of HST C/336 (blue), E6O6 (green), and Ji 2 s (red). The 
hnal column shows SExtractor segmentation maps of the pixels used for the photometry of each labeled galaxy component; the 
arbitrary color scale indicates component edges when two components are adjacent. Only pixel maps for the labeled objects 
are shown. Red (blue) circles (I'.'O diameter) indicate the centroid of the NB3420 emission (Lya emission). Photometry was 
performed individually on sub-arcsecond components associated with each LAE, and all components are labeled in the Veoe 
image. Postage stamps follow the conventional orientation, with north up and east to the left. The E—NB4670 color of each 
LAE is indicated below the object name, and objects marked by an asterisk (*) were found to have misidentified spectroscopic 
redshifts (see Section [57^ . 


tern for Eeoe, and the WFC3/IR DITHER-LINE pat¬ 
tern for J 125 and iLieo- In order to mitigate charge 
transfer efficiency (GTE) losses in our WFC3/UVIS 
Uf 336 exposures , we used the “post-fla sh” capability 
with FLASH =8 (|Biretta fc Baggettt I2013H . The final 3cr 
surface-brightness sensitivities and PSF FWHMs of the 
t^ 336 , Veoe, J 125 , and iLieo images are, respectively, 24.53, 
25.71, 25.79, and 25.56 mag arcsec-^ and 0"081, 0"092, 
0"178, and 0"186. The 3cr depths obtained in circular 
apertures with a diameter of 1.5 times the PSF FWHM 
are, respectively, 29.20, 30.22, 28.93, and 28.63 magni¬ 


tudes. 

4. DATA REDUCTION AND PHOTOMETRY 

Data reduction was performed on calibrated, flat- 
fielded, and CTE-corrected (in the case of WF C3/UVIS 
and ACS/WFC) images wi th DrizzlePac (jFrnchterl . 
120101: li^ekemoer et al.l . [20?)^ . The task Tweak Reg was 
used to align all exposures within each visit, and the 
AstroDrizzle pipeline was used to perform sky subtrac¬ 
tion, mask cosmic rays and bad pixels, and combine the 
exposures in the final, drizzled image. Final images were 
drizzled onto two scales: one optimized for the high- 
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NB3420 LyA I' F336W F606W F125W F160W UVJ Pixel Map 


Iae3038 

0.69 


Iae4070 

0.82 


Iae5200 

3.21 


laeSSIO 

0.68 



Figure 5. 5” X 5” postage stamp images of LAE photometric candidates (no spectroscopic confirmation) with NB3420 detections 
and imaging in all four HST filters. Objects are displayed and labeled as in Figure ID 


est resolution in each filter (for analysis of galaxy mor¬ 
phologies), and one where the pixel scale was consistent 
across all filters (for matched-aperture photometric anal¬ 
ysis). For the images drizzled to optimum resolution, 
the AstroDrizzle parameters pixscale and pixf rac are 
indicated in Table [2J For the images used in photo¬ 
metric analysis, all filters were drizzled to a pixel scale 
of 0.03 "/pixel, using a pixfrac value of 0.7 for C /336 
and 0.8 for Veoej >^ 125 , and F/ieo- The pixfrac values 
were chosen in order to achieve the optimum balance 
between the signal-to-noise ratio and visibility of low- 
surface-brightness features. Final drizzled images were 
registered to each other using the tasks TweakReg and 
TweakBack, achieving alignment between the Veos; J 125 , 
and Hiqo filters with an rms of 0"003 — O'.'OOb, and be¬ 
tween Veo6 and C /336 with an rms of O'.'Ol. The HST 
images were also aligned to the world coordinate system 
of the Keck/LRIS R—band image from M13 (the image 
to which all other Keck/LRIS images were registered) 
in order to map where the ground-based NB3420 (LyC) 
detections fell relative to emission in the HST images. 
After registration with TweakReg, residual astrometric 
distortions between the LRIS and HST images were cor¬ 
rected to a precision of 0"09 (less than half the size of an 
LRIS pixel) using the IRAF task CCMAP. Figures |3] — 
[S] display postage stamp images of the 16 galaxies with 
NB3420 detections and C/ 336 V 606 Ti 257 Li 60 imaging. 

For objects where imaging in all four filters 
(C/ 336 V 606 Ti 25 / 7 i 6 o) was available, the widest PSF was 
that of Hiqo (0"186). Accordingly, in order to perform 
matched-aperture photometry, we smoothed the higher- 
resolution HST images to match the PSF of the Hiqq 
image. When infrared data were unavailable because of 
the smaller field of view of the WFC3/IR instrument, the 
widest PSF was that of V/oe (0"092) and we smoothed 
the C /336 data to match this PSF. In order to perform the 
PSF-matching, we first created an empirical PSF from 
10—30 (depending on the filter) bright, isolated, and 
unsaturated stars using the IDL rout ine psf_extract 
from Star Finder (|Dioiaiti et al.L I2000D . These empiri¬ 


cal PSFs were then input into the IRAF routine PSF- 
MATCH, which outputs the convolution kernel and the 
PSF-matched image. The curves of growth of the stellar 
profiles in the iLieo and Veoe images agree with those of 
the PSF-matched images to <2% for the majority of the 
stellar profile, and agree within <10% at small radii (<3 
pixels). 

For objects with UqoqVqoqJ 125 Hiqo coverage, matched- 
aperture photometry wa s performed with SExtractor 
(jBertin fc Arnout^ . 119961) in dual-image mode, using the 
PSF-matched V 6 O 6 image to detect sources and define 
isophotes, and applying these isophotes (which can be 
examined using the SExtractor segmentation image) to 
the C/ 336 , Ti 25 , and iLi 6 o images. Because the Veoe im¬ 
age was already smoothed to the iLi6o PSF, no filter¬ 
ing was used, and the SExtractor detection threshold 
was set to 4.0cr. When photometry was extracted for 
galaxies without J 125 and Hiqq imaging, the unsmoothed 
Veoe image was filtered with a Gaussian kernel of cr=4 
pixels before source extraction. As most of the galax¬ 
ies in our z 2.85 sample have clumpy morphologies, 
SExtractor was run with maximum deblending (DE- 
BLEND_MINCONT = 0.0) in order to separate clumps 
within the galaxies for individual analysis. All galaxy 
clumps defined by SExtractor were examined visually in 
the segmentation image to guarantee that the subregions 
visible by eye were properly identified. For some galaxies, 
the detection threshold parameter was slightly modified 
(±lcr) to achieve the best isophote. There were also in¬ 
stances in which substructure was not visible in the Veoe 
image; in these cases, the relevant isophotes were defined 
in the C /336 or J 125 images where the substructure was 
visible. The final SExtractor segmentation maps used 
for photometry are displayed in the right-most column of 
Eigures[3] — [5l The background subtraction algorithm in 
SExtractor was set to LOCAL, which defines a square sky 
annulus around the object in question. However, in cases 
where multiple adjacent galaxy clumps were deblended, 
the annulus defined by SExtractor was often contami¬ 
nated by nearby sources to the extent that the resulting 
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background estimation was biased. In such cases, we cre¬ 
ated our own uncontaminated sky annulus around the 
object and estimated the background using the sigma- 
clipped mode, the procedure employed by SExtractor. 

In order to estimate photometric uncertainties for ob¬ 
jects i n the PSF-matched images, w e followed the meth¬ 
ods of iForster Schreiber et al.l (|2006ll and computed pho¬ 
tometric errors as a function of isophotal aperture size. 
First, we identified 1000 blank regions that avoided ob¬ 
jects and image edges in each of the PSF-matched t/ 336 , 
^ 606 ) >^125 j and Hieo images used for isophotal photom¬ 
etry. We then performed photometry on these blank re¬ 
gions with circular apertures of various sizes correspond¬ 
ing to the isophotal areas of our LBGs and LAEs. We 
defined the photometric uncertainty for a given object 
with isophotal area. A, as the standard deviation of the 
number of counts in the 1000 blank apertures of area A. 
The relationship between aperture size and background 
r ms in our images is quali t atively similar to that found 
in iForster Schreiber et al.l (|2006fi . From this analysis, 
we estimated 3a limiting magnitudes in apertures with 
a diameter of 1.5 times the PSF FWHM for each of the 
unsmoothed 1 / 336 , keoe, / 125 , and i/ieo images, and list 
them in Table[2] Photometric data for individual objects 
are listed in Table |31 

5. ESTIMATION OF PHOTOMETRIC REDSHIFTS 

One important challenge in isolating LyC emission 
from high-redshift galaxies is that the majority of these 
systems have complex morphologies. In our sample of 
z ^ 2.85 LBGs and LAEs, roughly 80% of the objects 
have complex morphology in the C/ 336 V 606 >/i 25 l/i 60 imag¬ 
ing, whether it be due to diffuse emission or multiple 
sources of nucleated emission. This fraction of objects 
with complex m orphology is similar to that found in 
iLaw et akl (|2007D for a sample of 66 z ~ 3 LBGs with 
rest-frame UV HST imaging 85%), and demonstrates 
that for most high-redshift galaxies, high-resolution im¬ 
ages reveal significant substructure. For the NB3420- 
detected galaxies, the NB3420 flux may be due to LyC 
emission from the high-redshift galaxy, or contamination 
from a lower redshift contaminant. Therefore, we must 
obtain photometric redshifts for individual subcompo¬ 
nents in order to identify possible contaminants. 

In this section, we begin by discussing the expected 
SED shapes of z ~ 2.85 galaxies and low-redshift con¬ 
taminants, and present the range of properties exhib¬ 
ited by galaxies in our sample. Next, we explain our 
procedures for fi t ting p hotometric redshifts with EAZY 
(jBrammer et al.l . 120081 1 to the objects in our HST sam¬ 
ple with imaging in all four of the C/ 336 , Veoe, J 125 , and 
i/i6o filters. Finally, we consider each NB3420-detected 
galaxy individually, discussing the photometric redshifts 
of each clump and the implications for the source of the 
galaxy’s NB3420 detection. 

5.1. Empirical Analysis of SEDs 

In order to describe the SED shapes of “typical” z ^ 
2.85 star-forming galaxies, we first present SEDs from 
the sample of LBGs and LAEs without LyC detections. 
This sample is described in more detail in Section 15.21 
Figure m shows t/ 336 V 6 O 6 / 125 /C 16 O photometry for several 
galaxies in the sample that span the range of typical SED 
shapes. As expected, the Lyman break is present in the 


SEDs of all galaxies without LyC detections. Further¬ 
more, every object in the LyC non-detection sample is 
undetected in C/336, which is consistent with the NB3420 
non-detections in M13. The presence of the Lyman break 
is one of the most important features for distinguishing 
low- and high-redshift galaxies in the sample of galaxies 
without LyC detections. It must be kept in mind, how¬ 
ever that the strength of the Lyman break may not be 
an effective way to identify the r edshi fts of LyC-emitting 
galaxies (see, e.g.. lCooke et al.l . I2014D . While the galax¬ 
ies in the LyC non-detection sample all exhibit a non¬ 
detection in C/336, there was a large range of V6O6 — /125 
colors in both the LAE and LBG samples, indicating the 
corresponding range in rest-frame UV slopes among the 
galaxies in our sample (see the upper panel of Figure [7]). 
In the LBG sample, the Veoe ~ /125 color ranged from 
—0.32 < Vqoq — J 125 < 1.03 (median Veoe ~ J 125 = 0-24), 
with most galaxies displaying red rest-frame UV slopes. 
In contrast, the LAE sample had —0.66 < Veoe — J 125 < 
1.7 (median Veoe — Ji 25 = —0.02), with most galaxies dis¬ 
playing blue rest-frame UV slopes. While on average the 
LAEs had bluer UV slopes than the LBGs, the LAE sam¬ 
ple also contained the galaxy with the reddest Veoe ~ J 125 
color (laelSfS, Zgpec = 2.847, which contains clumps with 
Veo6 — Ji 25 = 0-95 and V 606 ~/i 25 = 1.7). The Ji 25 ~Hieo 
colors of galaxies ranged from 0.20 < J 125 — i/ieo < 1-19 
(mean J125—i/ieo = 0.53), and did not differ significantly 
between the LBG and LAE samples (see the lower panel 
of Figure [7]). 

Additionally, there were a few objects in the LyC non¬ 
detection sample that demonstrated the typical SED 
shapes of low-redshift contaminants (see the bottom 
panel of Figure El). We found that the main features 
that help identify a low-redshift interloper are the pres¬ 
ence of flat C/336 ~ keoe and J125 — i/i60 colors, especially 
when accompanied by a large break between V6O6 and 
J 125 (see Figure ED. For contaminants, the break be¬ 
tween Veoe and J 125 corresponds to a Balmer break or 
4000A break between 0.5 < z < 2 . 2 , and has typical val¬ 
ues of 0.3 < Veoe — J 125 < 1-2. The only complication, 
however, is that young, dusty galaxies may also present 
nearly flat J 125 — i/ieo colors with red Veoe ~ J 125 col¬ 
ors. If C /336 “ heoe is flat (i.e., no Lyman break), then 
galaxies with flat J 125 — /Cieo colors and red Veoe ~ J 125 
colors can easily be identified as low-redshift contami¬ 
nants. However, if C /336 ~ heoe > 0, we must disentangle 
the degeneracy between contaminants and young, dusty 
galaxies. In our analysis, we label galaxies with such 
SEDs “ambiguous cases”, an example of which is shown 
in Figure El As the redshifts of galaxies with ambiguous 
SED shapes are uncertain, we do not include them in our 
analysis of z 2.85 galaxies without LyC detections. We 
also stress again that employing the C /336 “ heoe color to 
distinguish a contaminant from a z ~ 2.85 galaxy may 
be problematic when applied to galaxies with potential 
LyC emission, as the magnitude of the Ly man break for 
LyC galaxies is not well u nderstood fM13: iNestor et 
I2013I [Seidel et akl l2014fl . We keep this caveat in mind 
in Section 15.31 when discussing our targets with NB3420 
detections. Figure El demonstrates the relationships be¬ 
tween C /336 ~ ^ 06 j keoe ~ Jl 25 ^ and J 125 — Hieo colors 
for galaxies at z 2.85, galaxies with ambiguous SED 
shapes, and foreground contaminants. In the Veoe — J 125 , 
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Table 3 

LBG and LAE Photometry. 


ID*’' 

R.A.b 

Dec.*’ 

U336^ 

Vboe" 

J 125 " 

77i60" 

» d 

-^spec 

^phot 

^phot 

SED Type® 









PEGASE 

BPASS 


lael670a 

15:51:45.176 

19:10:15.261 

> 28.69 

29.39l°i® 

28.24121! 

> 28.64 

— 

1.19 

1.67 

c 

lael670b 

15:51:45.110 

19:10:15.225 

> 28.11 

28.25l°l° 

> 28.53 

> 28.18 

2.846 

2.28 

1.82 

r 

lae2292a 

15:51:47.633 

19:10:00.319 

> 29.01 

29.141°;!° 

28.971°°! 

28.6512;! 

2.851 

1.61 

1.44 

r 

lae2292b 

15:51:47.626 

19:10:01.034 

27.3ltg;?? 

07 70 + 0.10 

QQ 

26.7ll°;°! 

26.7612;!! 

— 

1.15 

1.17 

c 

lae2436a 

15:52:03.209 

19:12:51.261 

07 qq”I"0.32 
25 

26.98+!Jg? 

26.921°;“ 

on OO+ 0-08 
zu.oz_0 07 

2.04* 

2.29 

2.33 

c 

lae2436b 

15:52:03.239 

19:12:52.296 

26.131°;"^ 

25.161°;°! 

Z4.4tD_o 02 

04 OO+0.03 

z^.oo_o 03 

0.44* 

0.42 

1.65 

c 

lae2966a 

15:51:50.037 

19:10:42.064 

> 27.54 

26.821°;°! 

zu.zo_o 07 

26.2012;!2 

— 

1.09 

1.42 

c 

lae2966b 

15:51:49.975 

19:10:41.263 

> 27.14 

26.29l°;°6 

2fi 11 

ZD.li_o 08 

25.74l2;22 

2.841 

3.03 

2.90 

r 

lae3038a 

15:51:50.131 

19:09:02.336 

> 27.47 

28.021°;“ 

2fi 1 7+0.07 
ZD.l ^_0 07 

26.08l2;22 

— 

1.18 

3.10 

a 

lae3038b 

15:51:50.077 

19:09:02.483 

27.93to;2? 

28.201°;“ 

07 no+0-09 

z(.uo_o 09 

27.0812;!! 

— 

1.20 

1.21 

c 

lae4070a 

15:51:54.072 

19:10:26.789 

28.24t°;3l 

28.48l°l! 

28.3012“ 

07 7e:+0.19 
^(•(O_o.i6 

— 

1.25 

1.21 

a 

lae4070b 

15:51:54.061 

19:10:26.731 

07 qq-|-0.28 

z ( •00_Q 22 

27.5ll°;“ 

26.85l2;22 

26.8112;!! 

— 

1.24 

1.23 

c 

lae5200a 

15:52:01.125 

19:11:25.487 

OQ q'7-I-0.30 
_Q 23 

30.29l°;« 

> 29.74 

> 29.34 

— 

1.10 

0.99 

c 

lae5200b 

15:52:01.097 

19:11:25.994 

> 27.73 

27.90l°;!! 

26.69l2;2! 

26.5o12;2^ 

— 

1.23 

1.51 

a 

lae5200c 

15:52:01.111 

19:11:26.206 

> 27.84 

27.931°!° 

27.05l2;22 

27.4312;!! 

— 

1.11 

1.52 

c 

laefiSlOa 

15:52:07.004 

19:12:02.999 

> 28.68 

29.251°“ 

29.2612;!! 

> 28.90 

— 

2.03 

1.71 

c 

lae6510b 

15:52:06.976 

19:12:02.822 

> 27.71 

27.951°;“ 

27.6112;!! 

27.5812;!! 

— 

1.31 

1.35 

c 

laefiSlOc 

15:52:06.981 

19:12:03.350 

> 26.81 

26.591°!! 

0 *^ 70 + 0.06 
zo. / O_o 06 

25.58l2;2! 

— 

1.24 

1.38 

a 

lae6662a 

15:52:06.369 

19:10:42.590 

> 27.82 

07 “^1 ~l"0.09 

Z(.oi_o.o 8 

27.oi12;22 

26.64l2;22 

2.833 

2.82 

2.28 

r 

lae6662b 

15:52:06.350 

19:10:42.450 

> 27.49 

26.731°;°! 

26.68l2;22 

on 00 +O.O 8 

zu.zz_o 08 

2.833 

2.85 

2.80 

r 

lae6662c 

15:52:06.376 

19:10:43.706 

> 28.55 

28.841°;!° 

28.6312;!! 

28.6012;!| 

— 

1.26 

1.33 

c 

lae7180a 

15:52:04.669 

19:11:42.083 

> 27.62 

27.291°;!° 

Of\ 70+0-08 

zu. /z_o 08 

26.6712;!! 

— 

1.07 

1.54 

a 

lae7180b 

15:52:04.631 

19:11:42.088 

> 28.66 

29.191°;“ 

OQ qq+0.29 

2».50_o,23 

> 28.89 

— 

1.75 

1.74 

— 

lae7832a 

15:52:02.222 

19:10:48.807 

> 26.88 

26 . 081 °;°! 

26.08l2;2! 

25.6212;°° 

2.829 

2.77 

2.70 

r 

lae7832b 

15:52:02.199 

19:10:48.379 

> 27.91 

27.651°;“ 

27.6812 !! 

27.3712;!! 

2.829 

2.68 

1.86 

r 

lae7832c 

15:52:02.155 

19:10:48.614 

25.57l°;l? 

ot) 1 «+0.03 

zo.iu_Q Q3 

24.48l2;2! 

24.5i12;2| 

— 

0.72 

1.72 

c 

lae7890a 

15:52:01.957 

19:12:42.255 

26.6ll°;« 

2 fi 5O+0-06 

70 + 0.04 
zo. i o_0,04 

25.86l2;22 

— 

0.95 

1.38 

c 

M16a 

15:51:53.648 

19:09:29.392 

> 27.05 

26.15l°;°6 

Oc: 00+0.04 

zo.oz_o 04 

ot) 10+0-06 

zo.iy_o 05 

2.954 

3.92 

2.50 

r 

M16b 

15:51:53.619 

19:09:30.486 

> 27.89 

28.451°;“ 

27.2612;!! 

07 01 +0.19 
z(.zi_o 16 

— 

1.25 

1.36 

c 

MD5a 

15:51:45.206 

19:11:04.887 

> 27.21 

25.871°;°! 

25.681°;°° 

nc. 09 +O.O 6 

zo.oz_o 05 

3.143 

3.04 

2.88 

r 

MD5b 

15:51:45.226 

19:11:05.300 

> 27.05*5 

Oc: Qc+O.OS 

ZO.OO_Q Q4 

O’a ^7+0.06 

zo.<J(_o 05 

oc =0+0.08 

zo.oo_0.07 

3.143 

3.50 

1.94 

a 

MD12a 

15:51:51.887 

19:10:41.313 

> 28.70 

27.381°;°! 

26.92l2;2i 

on 51 +0.05 
ZD.Oi_o 05 

2.852 

3.17 

2.88 

r 

MD12b 

15:51:51.915 

19:10:41.281 

> 27.20 

2 fi C7+0-09 
ZD.Df_Q 08 

26.5912;!! 

26.1012“ 

2.852 

2.83 

2.56 

r 

MD12c 

15:51:51.866 

19:10:41.457 

> 27.15 

zo.uo_o,o4 

25.641°;°° 

9=91 +0.05 

zo.zi_Q Q 5 

2.852 

3.17 

2.91 

r 

MD12d 

15:51:51.879 

19:10:41.091 

> 27.56 

2fi 41 

ZD.4i_o 05 

26.3i12;2! 

25.97l2;2? 

2.852 

3.00 

2.94 

r 

MD12e 

15:51:51.880 

19:10:40.126 

27.291°;“ 

z(.zu_o 09 

26.62l2;22 

26.5712;!! 

— 

1.25 

1.25 

c 

MD12f 

15:51:52.028 

19:10:40.582 

> 27.81 

27 1 

ZMU_o.o7 

27.3512;!! 

07 0*^+®'^^ 
Z(.zo_Q 13 

2.852 

2.93 

2.74 

r 

MD12g 

15:51:52.036 

19:10:41.535 

27.261°;“ 

07 no“l"0.09 

z/.uo_o 09 

26.5o12;22 

on 1 5 + 0.10 

ZU.10_Q QQ 

— 

1.35 

1.26 

c 

MD12h 

15:51:52.035 

19:10:40.060 

> 28.86 

29.lll°;“ 

> 29.14 

> 28.79 

— 

1.95 

1.60 

— 

MD34a 

15:52:06.336 

19:12:48.673 

> 26.70 

25.491°;“ 

25.2112;°! 

24 n7+0-04 
Z4 .D^_q Q4 

2.852 

3.15 

2.87 

r 

MD34b 

15:52:06.307 

19:12:48.550 

> 26.83 

2 ^ 67"*”®'^^ 
ZO.Df_0.05 

24.841°;°° 

04 15 + 0.02 

Z4.10_q Q2 

2.852 

3.37 

2.66 

r 

MD34c 

15:52:06.334 

19:12:48.205 

> 27.60 

26.93lgg! 

26.151°;°° 

05 n7+0-04 
ZO.D^_q Q4 

2.852 

3.20 

2.65 

r 

MD34d 

15:52:06.314 

19:12:48.195 

> 29.54 

28.971°;“ 

28 . 061 °;°° 

07 4n+®-^^ 

Zf .4D_q Q5 

2.852 

2.71 

2.63 

r 

MD34e 

15:52:06.295 

19:12:48.191 

> 29.34 

28.941°;“ 

27.981°;°! 

07 qq+0.06 

4(.o5_q Q5 

2.852 

2.76 

2.65 

r 

MD34f 

15:52:06.318 

19:12:47.708 

> 27.27 

2fi QO+O-IO 
ZD.yz_o 09 

0*^ qq+0.05 

4o.50_o.o5 

05 45+0-05 

Z0.40_q Q5 

2.852 

2.12 

2.60 

r 

MD34g 

15:52:06.333 

19:12:47.380 

> 27.50 

z/.zo_o 10 

on 07 +O.O 6 

zu.z/_o 06 

on 90 +O.O 8 

ZU.ZO_Q Qg 

— 

1.21 

1.33 

c 

MD34h 

15:52:06.365 

19:12:47.446 

> 29.43 

29.901°;“ 

29.5412;“ 

29.3212;“ 

2.852 

1.73 

1.59 

a 


^ Objects are listed by their IDs from M13. The final letter of the object name indicates the sub-arcsecond component of the galaxy, according 
to the labels in Figures [3] [4] and [5] 

^ Object centroids calculated from Veoe- 

^ Isophotal magnitudes and photometric errors as a function of aperture size (see Section 0. Limits are 3cr. 

^ Spectroscopic redshifts listed were obtained via ground-based spectroscopy, and thus in general it is not possible to distinguish between 
individual sub-arcsecond components of galaxies within the spectrum. We list spectroscopic redshifts only when the SED fit indicates that the 
object is not a foreground contaminant. Spectroscopic redshifts are not listed for the following types of objects: LAEs without spectroscopic 
follow-up, contaminants identified via photometric redshifts, and objects where f /336 V606«^i25Lfi60 photometric data were insufficient to determine 
a photometric redshift. 

® This column indicates whether the SED of the obj ect i mplies a real 2 : ~ 2.85 galaxy (r) or a foreground contaminant (c). The letter (a) 
indicates an ambiguous SED shape, defined in Section l5.ll 

^ The original Lyo redshift associated with lae2436 from M13 was .2 = 2.832. Subseque nt rea nalysis of the available spectra near this object 
indicates that the 2 : = 2.832 emission was associated with a nearby Lya blob (see Section l5.3.1ll . 

s Although MD5b is formally undetected in U 33 Q at 3f7, emission at the location of MD5b is visible by eye in the t /336 imaging. This emission 
corresponds to a 2.25(7 detection in C /336 (^336 = 27.37^q'4q), which is consistent within errors of the detection in NB3420 (?T2ivs3420 = 
26.89l°;«). 
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MD7c 


Iae2158b 

(foreground 

contaminant) 



Figure 6 . Representative SED shapes found within our a 2.85 galaxy sample. The first three rows display z ~ 2.85 galaxies 
without NB3420 detections. From top to bottom, the objects shown are typical examples of SED shapes for objects with a red 
rest-frame UV slope (Veoe — •I 125 > 0.1), a flat rest-frame UV slope (—0.1 < Veoe — "I 125 < 0.1), a blue rest-frame UV slope 
(^606 ~ •I 125 < —0.1). The final row shows a typical example of the SED shape of a foreground contaminant. The left-most 
panel displays HST l/sse 1406^^125 composite color images, and indicates the sub-arcsecond components associated with each 
galaxy. The middle column shows the 17336 V606>li25l7i60 photometry for the relevant component (black circles with la error 
bars; downward-pointing arrows for Sct limits), the best-fit SED using the SMC-reddened BPASS models in EAZY (gray line), 
and the expected location of the photometric points based on the best-fit model (red circles). The right-most panel shows the 
redshift probability distribution for that component, with the soli d black curve indicati ng the probability distribution after a 
magnitude-based prior has been applied (see discussion of priors in iBrammer et ahl . 120081 1 and the dotted black curve indicating 
the probability distribution before applying the prior. Blue vertical lines indicate the best-fit photometric redshift (Zphot) and 
red lines indicate the observed spectroscopic redshift {zspec)- The top three objects demonstrate relatively narrow probability 
distributions that, while they do not have extremely high redshift precision, do encompass the true spectroscopic redshift. In 
general, the precision of the redshift probability distributions decrease with increasing photometric errors. The final object 
shows a clear example of a foreground contaminant with the photometric colors described in Section [5.II accordingly, the EAZY 
redshift probability distribution does not align with the spectroscopic redshift. For the case of lae2158, the actual LAE is 
lae2158a. 
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Jl25 ^160 


Figure 7. Histograms of Veoe — J 125 and J 125 — Hieo color for 
objects with Hsae 1406 ^ 125^160 photometry that do not have 
ambiguous SED shapes. LBGs are plotted in red, LAEs are 
plotted in blue, and contaminants are plotted in black. The 
top panel shows the wide range of Veoe — T 125 colors present in 
both the LAE and LEG samples, although on average LAEs 
have a bluer rest-frame UV slope than LBGs. The bottom 
panel demonstrates how important J 125 — Hieo color is in 
distinguishing typical z ~ 2.85 galaxies from contaminants, 
although there are some exceptions (e.g., objects with am¬ 
biguous SEDs; see open circles in Figure |8}. 

and J 125 — Hieo plot, z ^ 2.85 galaxies generally occupy 
a different region of color-color space from that of fore¬ 
ground contaminants, with the ambiguous cases (differ¬ 
entiated from contaminants by having U 33 Q — Veoe > 0 ) 
straddling both distributions. 

5.2. Modeling SEDs with EAZY 

While the empirical SED shapes suggest qualitative di¬ 
visions between low- and high-redshift galaxies, we can 
also obtain more systematic estimates of the redshifts 
of our targets. We used the photometric redshift code 
EAZY (|Brammer et al.l . [200811 to fit the C/ 336 , Veoe, J 125 , 
and Hiqq photometry acquired with El ST and estimate 
the redshifts of each sub-arcsecond component in the 
vicinity of our LAE and LEG targets. In order to choose 
the best input parameters for EAZY and learn how to 
interpret the output EAZY produces, we first ran the 
code on a test sample of z ^ 2.85 LBGs and LAEs 
in our El ST images with known spectroscopic redshifts, 
unambiguous SED shapes, and no NB3420 detections. 
Later, we supplemented this sample with z ~ 2.85 com¬ 
ponents of galaxies from the NB3420-detected sample for 
which there was no LyC detection (i.e., galaxies for which 
the NB3420 detection was proven to be associated with 
foreground contamination from another clump). This 



- 0.5 0.0 0.5 1.0 1.5 

J125 “ ^160 


Figure 8. U 336 — Veoe vs. J 125 — /C 160 and Veoe — 

J 125 vs. J 125 — /Cieo color-color plots of all objects with 
C /336 Veoe A 25 /C 160 photometry. LBGs with typical SED 
shapes are indicated by solid red circles, and LBGs with am¬ 
biguous SED shapes are indicated by open orange circles. 
LAEs with typical SED shapes are indicated by solid blue 
circles, and LAEs with ambiguous SED shapes are indicated 
by open light-blue circles. Foreground contaminants are in¬ 
dicated by black diamonds. Lower limits for all objects are 
indicated by upward-pointing arrows following the same color 
scheme. In the lower panel, the best candidate for true LyC 
emission (MD5b) is indicated by an orange open circle filled 
with green. 

test sample of galaxies should have a low rate of con¬ 
tamination by foreground interlopers, and thus help us 
evaluate whether or not EAZY can accurately identify 
galaxies known to be at z ~ 2.85 with the photometry 
provided. Additionally, we were able to analyze galax¬ 
ies in this sample that exhibited complex morphologies 
and evaluate the SED fits for each galaxy clump sepa¬ 
rately. In this way, we developed a procedure for differ¬ 
entiating between clumps that belonged to the spectro¬ 
scopically confirmed LBG or LAE and those that were 
lower redshift interlopers. As the galaxies in the test 
sample do not have NB3420 detections, we were able 
to make such a distinction without the complication of 
possible LyC emission, which might be associated with 
“non-standard” SEDs. Finally, we note that we did not 
include any components where the SED shape was am¬ 
biguous (see definition in Section [5Tj because it was not 
clear from the SED shape if that component was at high 
or low redshift. 

We varied several input parameters to EAZY in or¬ 
der to determine their optimal values. First, we experi¬ 
mented with fitting our data using different stellar pop¬ 
ulation synthesis models. EAZY defaults to PEGASE 
models (|Fioc fc Rocca-Volmerang^ Il997li . which span 
a range of star-formation histories, meta l licitie s, ages, 
and reddenings using the iCalzetti et al.l (|2000H atten¬ 
uation curve. To the five default PEGASE mod- 
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Figure 9. M16a demonstrates the typical shape of an “ambiguous” SED, one that can be fit with a wide range of redshifts 
that encompass both the spectroscopic redshift and lower redshifts typical of foreground contaminants. Results from EAZY 
are plotted in the three right hand panels, and both the SMC-reddened BRASS models (top row) and the PEGASE models 
that include nebular emission (bottom row) are shown to demonstrate the differences between models. In the left-most panel, 
the redshift is allowed to float during SED fitting, while in the middle panel the redshift is fixed to the spectroscopic redshift. 
The right-most panel shows the redshift probability distribution. Colors and symbols are as in Figure (6] The plot showing the 
PEGASE fit at fixed redshift demonstrates the limitations of the PEGASE models in terms of plotting young galaxies with dust 
attenuation; even with the redshift known, the PEGASE templates could not provide a satisfactory fit. Analysis of both the 
BRASS and PEGASE redshift probability distributions shows that for an SED of this shape, there is no way to narrow down 
the redshift to better than 0 < z < 4.5. 


els, EAZY adds additional templates representing a 
dusty starburst galaxy (described in iBrammer et al.l . 
1200811 and two templates for old, dusty galaxies (from 
EAZY_vl. l_lines/ and DUSTY/ in the online develop¬ 
ment versio rB). EAZY also includes a set of model tem¬ 
plates fro miBlanton fc RoweisI (1200^ (BR07), which are 
based on iBruzual &: Charlotl (|2003D models. Both the 
PEGASE and BR07 models provided good, qualitatively 
similar fits to most of our LBGs and LAEs, but failed 
to accurately represent galaxies with blue UV slopes 
(^606 — Ji 2 b < —0.1) and galaxies with SEDs younger 
than ^50 Myr (see Figure 0). Therefore, we experi¬ 
mented with additional stellar population models that 
might provide a better fit to the bluest and youngest 
galaxies in our sample. These included BPASS mod¬ 
els, which have a more accurate tre atment of Wolf-Rayet 
stars and massive stellar binaries (lEldridee fc StanwavL 
I2009D . and Starburst99 (S B99) models with up dated 
treatment of stellar rotation (jLeitherer et al.l . I2014D . Ad¬ 
ditional impetus for considering bluer templates is their 
increased emission in the LyC spectral region, which may 
provide a method of more accurately modeling galax¬ 
ies with LyC detections. In addition to experimenting 
with the choice of stellar population models, we also 
used two different extinction curves to redden the BPASS 
and SB99 models, for which only co nstant star-formation 
(CSF) histories are available. The iGalzetti et all (l2000f l 
attenuation curve has been traditionally used to model 
extinction in hig h-redshift st a r-forming galaxies, b ut re- 
cent work (e.g., iSiana etldl . 120091 : iReddv et al.l . I2006L 
[Mnl [Mil has shown that an SMC extinction curve 
may be more appropriate for galaxies with ages younger 
than 100 Myr. Therefore, we made two sets of CSF tem¬ 
plates with each of the BPASS and SB9 9 models, one 
template set reddened exclusively with the lCalzetti et al.l 

® https://github.com/gbrammer/eazy-photoz/blob/master/ 

templates/ 


(l200r)tl attenuation cury e, and another reddened with 
the lCalzetti et al.l (12000( 1 attenuation curve for templates 
with older galactic ages a nd the SMC extinction curve 
from iCordon et al.l (1200,8(1 for templates with ages less 
than 100 Myr. 

Based on the analysis of objects without LyC de¬ 
tections, we verified that, for the majority of galaxies, 
EAZY estimates the correct redshift of the galaxy within 
roughly ±0.5 of the spectroscopic redshift. Thus, we 
can successfully use photometric redshifts to determine 
if most galaxy clumps are at z ~ 2.85 or are low-redshift 
contaminants, keeping in mind that we have fairly coarse 
redshift precision. As EAZY provides several estimators 
for the photometric redshift (e.g., ZpeaJEl), 

we investigated each of them while varying the input pa¬ 
rameters to the program and determined that Zpeafe pro¬ 
vided the best estimate of the spectroscopic redshift. All 
photometric redshifts we quote use the Zpeak estimator. 
Of the range of input parameters tested on the sample 
of objects without LyC detections, we determined that 
the PEGASE models result in photometric redshifts that 
most closely match the spectroscopic redshifts. These 
models had both the smallest systematic offset and low¬ 
est standard deviation of all the variations of input pa¬ 
rameters we tried, yielding {zphot - ,2spec)/(l ± -^spec) — 
—0.03 ± 0.07. We note that while the PEGASE models 
work best for the sample as a whole, the BPASS mod¬ 
els provide the best fits to galaxies with blue rest-frame 
UV slopes (Veoe — Ji 2 b < —0.1). Also, as the EAZY 
PEGASE models only include 5 templates, they did not 

® Redshift where the likelihood is maximized after applying the 
magnitude-based prior. 

Redshift marginalized over the posterior probability distribu¬ 
tion. 

Hybrid between Zp and Zm 2 to address the pathological case 
where there are two widely-separated peaks in the probability dis¬ 
tribution that have similar integrated probabilities. 
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accurately fit galaxies with younger SEDs. Better fits 
were achieved for young galaxies by implementing a fine 
grid of young BRASS templates, specifically by includ¬ 
ing SMC-reddened models with ages of 1 Myr, 5 Myr, 

10 Myr, 30 Myr, and 50 Myr. Throughout our analysis 
of the galaxies with potential LyC detections (Section 
15.31) , we employ both the PEGASE models and the SMC- 
reddened BRASS models with additional young galaxy 
templates in order to fully examine the likely photomet¬ 
ric redshifts for each galaxy. All figures in this paper 
displaying output from EAZY show fits using the SMC- 
reddened BRASS models, unless otherwise indicated. 

5.3. Results of Photometric Redshift Fits for Potential 
LyC Emitters 

In total, we have observations in all four 
£^ 336 ^ 606 ^^ 125^^160 filters for 16 galaxies with poten¬ 
tial LyC detections. These 16 galaxies include 4 LBGs 
and 8 LAEs with spectroscopic confirmation, as well 
as 4 LAE photometric candidates. In this section, we 
discuss the results of our analysis of these potential LyC 
emitters. We describe the sources of contamination for 

11 targets with obvious contaminants, present 4 galaxies 
with ambiguous SEDs that may or may not be at high 
redshift, and argue for MD5 as a true LyC emitter. We 
also note 3 objects for which the HST imaging revealed 
that the spectroscopic redshift was incorrectly assigned 
in M13. This analysis is contingent upon the precise 
alignment between the HST and LRIS NB3420 imaging, 
as described in Section SI 

5.3.1. Foreground Contaminants in the LyC Sample 

Seven objects with Keck/LRIS spectroscopic redshifts 
(2 LBGs: MD12, M16; 5 LAEs: lael670, lae2292, 
lae2966, lae6662, and lae7832) and one object in the pho¬ 
tometric LAE sample {lae3038) had contaminants that 
stood out plainly with the combination of high-resolution 
HST imaging and SED fitting. Each of these objects was 
resolved into several clumps in the HST imaging. In all 
cases, at least one of the clumps had an SED fit corre¬ 
sponding to the redshift of the Keck/LRIS spectrum (or 
consistent with the spike redshift z = 2.85, in the case of 
lae3038), while the clump associated with the LyG emis¬ 
sion had the unambiguous SED of a 0.5 < z < 2.2 con¬ 
taminant (flat Ussq — Vqoq and Ji 25 ~Hi 60 j red V 606 ~>^i 25 ; 
similar to the example contaminant shown in the bottom 
panel of Figure |6]). For these eight objects, the NB3420 
emission had a fairly large offset from the original galaxy 
coordinates; the offsets between the NB3420 detection 
and Lya emission ranged from 0V57 to 1'.'15 with a me¬ 
dian value of 0'.'65 (5.0 kpc at z = 2.85), and offsets 
between the NB3420 emission and the LRIS R—band 
emission ranged from 0'.'12 to I'.'26 with a median value 
of 0'.'72 (5.7 kpc at z = 2.85). The morphology of these 
objects in the HST images supports the evidence from 
the SED fits that these objects are not physically associ¬ 
ated; the clumps are distinct, with no evidence for diffuse 
emission between them that might indicate interactions 
between galaxies at the same redshift. 

A third LBG in our sample, MD34, has an NB3420 
detection coincident with two clumps in the HST imag¬ 
ing, only one of which is from a foreground contaminant. 
Figure |3] shows the eight distinct components associ¬ 
ated with MD34 (MD34a though MD34h). Gomponents 


MD34a through MD34f all have SED shapes that place 
them at the spectroscopic redshift of MD34 (z = 2.85), 
and MD34g and MD34h have SED shapes indicative 
of foreground contaminants. The NB3420 emission for 
MD34 is coincident with two clumps in the HST imag¬ 
ing: MD34f (z = 2.85) and MD34g (z < 2.5), separated 
on the sky by 0'.'31. Both MD34f and MD34g exhibit 
emission in the C /336 filter, although both detections are 
less than 3 (t. It is unclear whether the NB3420 emission 
associated with MD34 is due solely to non-ionizing UV 
radiation from the foreground contaminant MD34g, or 
also in part to LyG emission from the z = 2.85 compo¬ 
nent MD34f. As the high-resolution U 33 Q image does not 
exclusively probe the LyC spectral region at z = 2.85, 
and the foreground contaminant MD34g is too close to 
MD34f to distinguish in the seeing-limited NB3420 im¬ 
age, we cannot confirm MD34 as having a robust LyC 
detection. Because the interpretation of this case is am¬ 
biguous, we do not include MD34 in the final sample of 
LyC emitters. 

For an additional object, lae2436, the prese nce of an 
extended Lya blob (|Steidel et al.l . 1200(1 1201 Ih near the 
position of the LAE made the original analysis of the 
LRIS imaging and spectra difficult. The deep, high- 
resolution HST imaging helped us clarify the interpreta¬ 
tion of this object and determine that the NB3420 detec¬ 
tion is associated with a foreground contaminant. Figure 
[iQ] shows the complex morphology of lae2436 in both 
the LRIS R—band and the continuum-subtracted Lya 
images. The Lya blob extends over more than 5", and 
in the LRIS R—band image two bright galaxies (X and 
Y) appear to be in the vicinity of the Lya blob, with 
diffuse emission between them. As shown in Figure 1101 
both Galaxy X and Galaxy Y exhibit Lya emission in 
the LRIS spectrum (L3 and LI, respectively) and were 
identified as LAEs in M13. Galaxy X was originally iden¬ 
tified as the LAE lae2436 (z = 2.832), and this galaxy 
is clearly detected in the NB3420 image. There are no 
other NB3420 detections nearby. Finally, there is a Lya 
emission line (L2) in the LRIS spectrum coincident with 
the diffuse emission between Galaxies X and Y in the 
LRIS R—band image. This diffuse emission was orig¬ 
inally attributed to the presence of the Lya blob and 
thought to be associated with Galaxy X, as L2 and L3 
have nearly identical wavelengths. 

In the HST imaging, much additional substructure is 
revealed in the vicinity of the Lya blob. First, it becomes 
clear that the diffuse emission in the LRIS R—band im¬ 
age is due to several unresolved, faint galaxies in close 
proximity. One of these galaxies (Z) is roughly coincident 
with the L2 Lya emission line in the LRIS spectrum, and 
the SED fit to this galaxy using HST photometry places 
it near z = 2.832, the redshift of the emission lineQ 
None of the SEDs for other galaxies in the vicinity of 
the L2 emission line demonstrate the typical features of 
z ^ 2.85 galaxies. As for object X, originally identified as 
lae2436, it can be distinguished as two separate galaxies 
in the HST image (a and b), both of which are detected 
in NB3420. Subsequent reanalysis of the available spec- 

We note that Galaxy Z has a unique and extremely red SED, 
with Veo6 — J 125 = 2.08 and J125 — Hieo = 1.18. This galaxy is 
a sub-millimeter source that will be further described in Steidel et 
al., in prep. 
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Figure 10. HST and LRIS Imaging for lae2436 is dis¬ 
played (9" 5 X 9" 5), along with the LRIS spectrum prob¬ 
ing the rest-frame UV at z ^ 2.85. Along the right-hand 
column of the figure, the morphology of the Lya-blob in the 
vicinity of lae2436 is shown in LRIS V (non-ionizing UV con¬ 
tinuum), LRIS continuum-subtracted NB4670—U (indicating 
Lya emission), and LRIS NB3420 (LyC emission). To the left 
of the LRIS V image, the higher resolution HST Veoe image 
is shown. Blue contours indicate the location of the Lya 
emission based on the NB4670—U image. The red rectangle 
in the HST Veoe image shows the location of the l”2 slit, 
and the LRIS 2D spectrum below is aligned to match in the 
orientation and spatial scale of the imaging. The zoomed-in 
portion of the 2D spectrum shows the Lycr emission (LI, L2, 
L3), and the full 2D spectrum is displayed on the left. In the 
original LRIS imaging, only two galaxies were visible near the 
Lya-blob (X and Y). These galaxies were identified as LAEs 
based on their 1/—NB4670 colors and the Lya emission in 
the LRIS spectrum. As discussed in Section f5. 3. II Galaxy Z 
was first identified in the HST imaging and is likely associ¬ 
ated with the Lya emission in the LRIS spectrum. Galaxy 
X, meanwhile, can be distinguished as two segments in the 
HST image (a and b) that were both identified spectroscop¬ 
ically as foreground contaminants through reexamination of 
the available LRIS spectra. The [Oil] {z = 0.44) emission 
from lae2436b is visible in the top of the full 2D spectrum. 
We conclude that there is no NB3420 emission associated with 
lae2436. 

tra near the Lya blob allowed us to confirm spectroscopic 
redshifts for both of these galaxies {z = 2.04 for lae2436a, 
z = 0.44 for lae2436b), and the SED fits to both objects 
are consistent with their spectroscopic redshifts. In par¬ 
ticular, we draw attention to the spectrum of lae2436b, 
visible in Eigure [10] on the right-hand side of the LRIS 
spectrum, coincident with L3. In the zoomed-out ver¬ 
sion of the spectrum, a spurious emission line is visible 
at 5375 A which we have identified as [Oil] at z = 0.44. 
The Lya emission line originally identified for this object 


(L3) must be due to extended Lya emission from the 
Lya-blob. In conclusion, as both lae2436a and lae2436b 
are at z < 2.82, the NB3420 filter does not probe LyC 
emission for these galaxies and thus the NB3420 detec¬ 
tions associated with both galaxies are foreground con¬ 
tamination. 

Einally, we found that one object in the spectroscopic 
LAE sample, lae7890, was misidentified as an LAE. This 
object was presented in the Appendix of M13 as a faint 
LAE with a borderline color excess (U—NB4670 = 0.70) 
and a possible LyC detection, and was not analyzed with 
the main LAE sample for which a NB4670 magnitude 
limit of 7714670 < 26 was imposed. A marginal emission 
line had been identified for this object, placing it at the 
spike redshift of z = 2.85, but the SED (similar to that 
of the contaminant in Eigure |6|) indicates unequivocally 
that lae7890 is at low redshift, in the range of 1 < z < 2. 
We considered [Oil] as a possible source of the emission 
line, but that would place lae7890 at z = 0.26, which is 
also inconsistent with the observed SED. Thus we con¬ 
clude that this faint emission line is either spurious, or 
is possibly consistent with [CIII] I907A / CIII] I909A 
emission from a galaxy at z = 1.45. In either case, the 
NB3420 emission is not LyC. 

In summary, II/I6 candidate LyC emitters in our sam¬ 
ple show obvious signs of foreground contamination at 
the position of the NB3420 detection. Eor 9 objects, 
there is a bona fide z = 2.85 galaxy at the position of the 
LAE, with an additional foreground galaxy offset from 
the LAE and associated with the NB3420 emission. Eor 
2 objects, z = 2.85 was erroneously assigned to the can¬ 
didate LAE, and, again, the NB3420 emission is actually 
non-ionizing UV flux from a low-redshift contaminant. 
In all of these cases, the available evidence suggests that 
we are not observing LyC emission at z = 2.85. 

5.3.2. Ambiguous Cases 

Eor the four photometric LAE candidates with NB3420 
detections, the SED shape was ambiguous and there were 
no spectroscopic redshifts available to confirm that the 
objects are indeed at z ~ 2.85. These two factors make it 
impossible to confidently claim a LyC detection for any 
of these objects. We discuss the photometry and SED 
fits for these objects in detail in Appendix A, and sum¬ 
marize the results here. One object, lae4070, has similar 
J 125 — H 16 O and Veoe —Ti 25 colors to many z ~ 2.85 galax¬ 
ies in our sample and is the most promising photometric 
LAE candidate for true LyC emission. The other three 
objects {lae5200, lae6510, and lae7180) disp lay the am¬ 
biguous SED shape described in Section [ST] which may 
describe galaxies at many redshifts. As we cannot unam¬ 
biguously determine whether or not the four photomet¬ 
ric LAE candidates discussed in this section are truly at 
z ~ 2.85, we adopt a conservative approach and do not 
count the NB3420 detections for these objects as secure 
signatures of leaking LyC radiation. 

5.3.3. LyC Emission from MD5 

The best candidate for true LyC emission is the LBG 
MD5, which has a spectroscopic redshift of z = 3.14 con¬ 
firmed by spectra from both LRIS and the MOSFIRE 
near-IR multi-object spectrograph (see Figure (TTI) . The 
LRIS spectrum was taken in May 2011 (MI3) and shows 


















A High-Resolution HST Study of Apparent LyC Leakers at z ^ 3 


15 


MD5 


1 

) 

——— 

a 




5000 5100 5200 5300 5400 5500 

Wavelength [A] 


MD5a 



MD5b 



i Fixed z (z= 

s i- 

7 i- 

3.14) 

10^ 10^ 1 
wavelength [A] 

: Fixed z (z= 

s i- 

7 i- 

,1_J 

3.14) 




Figure 11. Imaging, spectra, photometry, and SED fits for MD5. The top left panel shows the 17336 ^ 606^25 color-composite 
HST image for MD5, indicating the two sub-arcsecond components (MD5a and MD5b). The locations of the MOSFIRE (yellow) 
and LRIS (magenta) slits are overlaid on the image. While both MD5a and MD5b fall within the MOSFIRE and LRIS slits, it is 
not possible to distinguish between the clumps in the spectra due to the slit orientation and seeing (076 for MOSFIRE, 077-l70 
for LRIS). The ID MOSFIRE and LRIS spectra of MD5 are displayed to the right of the image, showing superimposed emission 
from both sub-arcsecond components. We detect the [OIII] doublet (4959A, 5007A; z = 3.1426) in the MOSFIRE A—band 
spectrum, and Lya emission {z = 3.147) along with interstellar absorption lines (1303A(OI -I- Sill), 1334A(CII); z = 3.139) in 
the LRIS spectrum. The interstellar absorption lines are blueshifted relative to the Lya and [OIII] emission, indicative of an 
outflow. Below the spectra and imaging are shown the EAZY output for MD5a and MD5b. In the left-most panel, the redshift 
is allowed to float during SED htting, while in the middle panel the redshift is Hxed to z = 3.14, the redshift indicated by the 
spectra. The right-most panel shows the redshift probability distribution. Colors and symbols are as in Figure [6l MD5a has 
the typical SED shape associated with z ~ 3 galaxies in our sample. MD5b, which is associated with the NB3420 detection, 
exhibits the ambiguous SED shape discussed in Section [5. II As there is no evidence of foreground contamination in either the 
MOSFIRE or LRIS spectra, we propose that both of these clumps are at the spectroscopic redshift 3.14 and that the NB3420 
emission associated with MD5b is true LyC emission. We note that the LRIS and MOSFIRE spectra rule out strong emission 
lines in several redshift ranges, including Ha emission between 1.97 < z < 2.65, Lya emission at z > 1.5, [Oil] emission at 
z < 0.9, and [OIII] emission at z < 0.4. 


Lya emission (z = 3.147), along with multiple absorption 
features (CII, Sill, and 01; z = 3.139). The AT—band 
MOSFIRE spectrum was acquired in June 2012, and con¬ 
tains detections of both [OIII] emission lines (4959A, 
5007A; z = 3.1426). Unfortunately, we were unable to 
measure the z = 3.14 H/3 emission line in this spectrum 
because it falls on a sky line. In the HST imaging, MD5 
is composed of two clumps separated by 0758 along a di¬ 
rection 32 degrees East of North (MD5a and MD5b; see 
Figure ITT]) : only MD5b is coincident with the NB3420 
emission. The orientation of the MOSFIRE and LRIS 
observations (slit PAs of 274 and III degrees, respec¬ 
tively; see Figure HH) were such that both clumps fell 
within the spectroscopic slits. Although it is not possi¬ 
ble to distinguish between the clumps in the spectra due 
to slit orientation and seeing (076 for MOSFIRE, 077-170 
for LRIS), neither spectrum shows evidence for spurious 
emission or absorption features that would indicate the 


presence of a low-redshift interloper _ 

Ideally, as in the case of lonl in iVanzella et al.l (|20I2f) . 
a candidate LyC-emitting galaxy would have a simple, 
compact morphology. In such a case, the probability 
of a foreground interloper is negligible. However, it 
has been shown that high-redshift gal axies typica l ly ex¬ 
hibit clumpy mo rphologies (see, e.g., iLaw et ahl . I2007t 
iLotz et all . I2006D . Of the 35 LBGs in our z > 2.82 sam¬ 
ple with imaging in Veoej only 20% have simple, com¬ 
pact morphologies. Close inspection reveals clumpy mor¬ 
phologies in all other cases, and several LBGs are com¬ 
prised of clumps with offsets significantly greater than 
the 0758 offset of the MD5 clumps. For example, LBGs 
M23, MD9, and CI3 are comprised of clumps suggested 
by SED fits to be at the LBG redshift with offsets of 0791, 
I700, and 1719, respectively. As the majority of LBGs 
consist of several associated clumps at the same redshift, 
the fact that MD5 displays multiple clumps does not nec¬ 
essarily indicate that MD5b (the clump associated with 
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the NB3420 detection) is a foreground contaminant. 

However, the presence of the second clump opens the 
possibility that one of these clumps is a low-redshift in¬ 
terloper, and thus we examine the SED fits to MD5a and 
MD5b for evidence of foreground contamination (see Fig¬ 
ure ni). Both clumps are nearly identical in Vqqq mag¬ 
nitude, and are within 0 . 1 - 0.2 magnitudes in J125 and 
iJi 6 o- MD5a has the SED shape typical of most z ^ 3 
galaxies in our sample and is almost certainly at the spec¬ 
troscopic redshift z = 3.14. MD5b has the ambiguous 
SED shape described in Section 15.11 As the spectro¬ 
scopic redshift of MD5 is higher than most galaxies in 
our sample, the Hieo filter actually falls right at the lo¬ 
cation of the Balmer break, rather than redwards of the 
break. Thus, the i?i 6 o filter is partially contaminated by 
flux bluewards of the Balmer break, and J 125 — i?i 60 no 
longer probes the full strength of the break. Therefore, 
the fact that MD5b has a smaller J 125 — Hieo color than 
the lower-redshift LBGs in our sample does not neces¬ 
sarily indicate a young age or low-redshift contaminant. 

Figure [ 8 ] shows MD5b plotted with respect to other 
galaxies in our sample in Veoe ~ <^125 vs. J 125 — l?i 60 - If 
MD5b is truly at z = 3.14, it has a very small Balmer 
break (J125 —H^ieo = 0.04±0.13), a red UV slope (Veoe — 
•^125 = 0.28 ±0.10, equal to the mean of the LBG sample), 
and a C /336 detection brighter than expected by standard 
stellar population synthesis models (unsurprising for a 
LyC emitter). If MD5b is a foreground interloper, EAZY 
estimates that the most likely redshifts for the interloper 
are 0 < z < 0.5 or 2 < z < 2.5 using PEGASE models, 
and between 1 < z < 2.8 using reddened BPASS models. 
In both cases, z ~ 2.3 —2.4 is the most likely contaminant 
redshift. 

The strongest evidence that MD5b is indeed at z = 
3.14 is the lack of spurious emission or absorption lines 
in both the LRIS and MOSFIRE spectra. Because these 
spectra span different wavelength ranges (3100—7000 
A for LRIS and 1.95—2.4 ^m for MOSFIRE), we can 
rule out strong emission lines for several redshift ranges. 
These include Ha emission between 1.97 < z < 2.65, 
Lya emission at z > 1.5, [Oil] emission at z < 0.9, and 
[OIII] emission at z < 0.4. We note that the lack of 
spurious emission lines in the AT-band MOSFIRE spec¬ 
trum rules out Ha emission right in the redshift range 
predicted by EAZY to be the most probable redshift of a 
contaminant (2 < z < 2.6). It might be possible to con¬ 
firm the redshift of MD5b using MOSFIRE observations 
with the slit oriented along the axis connecting MD5a 
and MD5b. These observations would maximize the dis¬ 
tance between the two clumps (separated by 0758) and, 
if taken under conditions of good seeing, potentially dis¬ 
tinguish emission from each clump individually. 

6 . PROPERTIES OF THE LYMAN-CONTINUUM EMITTER 

MD5 

One of the main goals of this work is to investigate the 
multiwavelength properties of galaxies with and without 
LyC emission, in order to better understand the mecha¬ 
nism of LyC photon escape from galaxies. We are also 
interested in investigating any systematic differences be¬ 
tween galaxies with and without LyC emission, for such 
differences may facilitate the search for LyC-emitting 
galaxies both during and after the epoch of reioniza¬ 
tion. As our HST data has left us with only one robust 


candidate for LyC emission (MD5), we here discuss the 
morphological properties and best-fit stellar population 
of this object with respect to the properties of typical 
LBGs. 

6 . 1 . Morphology of MD5 

Figure [3] displays imaging of MD5 in all available 
bands. Morphologically, MD5 is composed of two clumps 
(MD5a, MD5b) separated by 0758 (4.4 kpc at z = 3.14). 
The Veoe magnitudes for each clump (representing the 
non-ionizing UV continuum) are = 25.87 ± 0.04 

and = 25.85lQ g 4 . I'^ Section 15.3.31 we present 

arguments for why both of these clumps are likely at the 
spectroscopic redshift of z = 3.14 and why foreground 
contamination is unlikely. The clumpy morphology of 
MD5 is similar to that of many other LBG systems, 
which commonly exhibit significant substructure. MD5b, 
which is more diffuse and lower in surface brightness than 
MD5a, is the clump associated with the NB3420 detec¬ 
tion (i.e., the LyG emission). 

Because of the particularly high redshift of MD5, we 
are able to directly map the LyC emission in the high- 
resolution HST 1/336 image. While most galaxies in our 
sample had redshifts of roughly z ~ 2.85 and the t /336 fil¬ 
ter was partially contaminated by non-ionizing flux red- 
wards of the Lyman limit, MD5 is at high enough redshift 
(z = 3.14) that the t /336 filter probes the LyC spectral 
region without any contamination. While MD5b is for¬ 
mally undetected in t /336 at 3cr, emission at the location 
of MD5b is visible by eye in the t /336 imaging. MD5b 
has a 2.25cr detection in t /336 (^335 = 27.37tQ 4 Q), which 
is consistent within errors of the detection in NB3420 
{mNB342o = 26 . 891 q ; 34)0 

In light of models in which LyC emission may escape 
anisotropically from g alaxies (e.g., I Gnedin et al.l . I2008t 
iZackrisson et al.l . 1201311 . we examined the offset between 
the centroid of the t /336 and V 6 O 6 emission for MD5b in 
order to determine if there was a significant offset be¬ 
tween the ionizing and non-ionizing UV emission. We 
measured a value for this offset of Ajjv = 0708. As there 
are no additional z ~ 3 galaxies in our sample with t /336 
detections to compare to, we examined the distribution 
of t /336 ~ keoe offsets for the foreground contaminants in 
our sample. In this way, we measured A[/y for objects 
where t /336 and Veoe are both probing the non-ionizing 
continuum, and thus should not demonstrate significant 
offsets. For the contaminant sample, we found a roughly 
flat distribution of offsets between 070 and 0712 with a 
mean offset of 07065 and a standard deviation of 07031. 
As MD5b has an offset consistent with the mean of the 
contaminant distribution, we conclude that its measured 
offset is not significant. This lack of significant offset im¬ 
plies either that LyC emission is escaping isotropically 
from MD5b, or that, if LyC emission escapes only from 
cleared holes in the ISM, MD5b must be oriented such 
that the opening is along our line of sight. 

6 . 2 . ISM Kinematics of MD5 

The spectral features of MD5 shed light on the kine¬ 
matics of its ISM. The [OIII] A5007 nebular emission line 

We note that the ±336 filter is wider than NB3420, and thus 
the fainter f/336 magnitude may be due to increased IGM attenu¬ 
ation within the bluer half of the Usse filter. 
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observed in the MOSFIRE spectrum of MD5, which in¬ 
dicates the systemic redshift, places MD5 at z = 3.1426. 
This line has an intrinsic width of cr„ = 37 km s“^, typ¬ 
ical of LAE s, but half that of t ypical LBGs (Trainor et 
ah, in oreo: IPettini et l2nnih . The low-ionization in¬ 
terstellar absorption features observed in the LRIS spec¬ 
trum of MD5 (01 -I- Sill A1303, CII A1334) are con¬ 
sistent with a redshift of z = 3.139, blueshifted relative 
to the [OIII] emission. The magnitude of this blueshift 
corresponds to a velocity offset of Avjs ^ 280 km s“^, 
higher than t h e med ian Avjs for LBGs (150 km s“^; 
iShaolev et al.l . 1200311 . but fairly uncertain due to the 
low signal-to-noise of the LRIS spectrum. Addition¬ 
ally, the redshift derived from the centroid of the Lya 
emission line (z = 3.147) corresponds to a velocity off¬ 
set of AvLya ~ 300 km s“^, which is typical of LBGs 
{Av^y^ = 360 km s~^: IShaolev et all . I2003fl . Thus, in 
terms of its kinematics, MD5 does not stand out sig¬ 
nificantly with respect to the full population of LBGs, 
although the blueshift of its interstellar absorption lines 
is higher than average. We note that with the spatial 
resolution of our spectra, we cannot resolve the two com¬ 
ponents MD5a and MD5b separately and evaluate their 
individual kinematics. Finally, we note that our mea¬ 
sured value of AvLya is in consistent with predictio ns for 
LyG-leaking galaxies from iVerhamme et al.l (j2014l) , who 
find small offsets between Lya emission and the systemic 
redshift ( AvLya < 150 km s“^) in models of galaxies with 
leaking LyC radiation. 

6.3. Stellar Populations of MD5 

In addition to providing high-resolution imaging of 
MD5, the multiwavelength HST data also enabled us 
to fit the photometry for both sub-arcsecond compo¬ 
nents (MD5a and MD5b) with stellar population syn¬ 
thesis (SPS) models. While EAZY is a powerful tool for 
estimating photometric redshifts, the program does not 
directly provide stellar population information. In order 
to examine the stellar populations of the galaxies in our 
sample , we employed the stellar population fitting code 
FAST (|Kriek et al.Ll200^ . which accommodates a differ¬ 
ent set of stellar population synthesis models and dust 
extinction prescriptions from those of EAZY. To model 
the photometry, we fixed the r edshift to the spectr o scopi - 
cally measured value and used iBrnzual fc Charlotl (j20n3[ i 
models ranging in age from 50 Myr to 2 Gyr with delayed- 
T star -formation histories (SFR oc te~*^'^) and a lChabri^ 
(|2003li initial mass function. We chose delayed-r star- 
formation histories because of their flexibility in accom¬ 
modating both rising and falling star-formation histo¬ 
ries, and we note that we found no significant qualita¬ 
tive differences in the derived stellar population param¬ 
eters when experimenting with constant or rising star- 
formation histories. The lower age limit of 50 Myr is 
ado pted to reflec t the L BG dynamical timescale, follow¬ 
ing iReddx^eF^ (I2012D . However, this age limit is con¬ 
servative and may be larger than necessary given the 
small sizes of some galax y subcomponen t s. Fo r dust ex¬ 
tinction, we employed the iCalzetti et al.l (j200nfl attenua¬ 
tion curve. We adopted solar metallicity for the models 
because, with only four photometric points, we did not 
have enough data to constrain metallicity. By perform¬ 
ing tests with several values for fixed metallicity and with 


Rest-Frame Wavelength [Ltml 
0.12 0.24 0.48 li 2A 



Figure 12. Best-fit IBruzual fc Charlotl (120031! stellar pop¬ 
ulation synthesis models as computed by FAST (gray line) 
to HST 1 / 336 F 606 Ji 2577 i 60 photometry for MD5a and MD5b 
(black circles), fit at the spectroscopic redshift of z = 3.14. 
One-sigma photometric uncertainties are smaller than the 
data points. For MD5b, the 2.25a detection in Usse (probing 
the LyC spectral region at z = 3.14 with no contamination 
redwards of the Lyman limit) is plotted in green, with the for¬ 
mal 3a limit indicated by the black arrow. We note that while 
the 1/336 data point for MD5b does not a gree with the model 
predic tion, this is to be expected from IBruzual fc Charlotl 
(I 2 OO 3 I! models, which do not have significant emission blue- 
wards of 912 A. Parameters of the fits are listed in Table [d] 

metallicity as a free parameter, we verified that metallic¬ 
ity does not have a significant effect on the final fit. We 
note that even when metallicity was allowed to float, the 
best-ht model for MD5b had solar metallicity. 

We display the best fit models and photometry for 
MD5a and MD5b in Figure [T^ We note that while the 
C ^336 data point for the LyG-emitter MD5b does not agree 
with the model predictio n, this is to be expected from 
IBruzual fc Gharlotl l|2003ll models, which do not have sig¬ 
nificant emission bluewards of 912A. We also report the 
best-ht values and 68% conhdence intervals for r, stel¬ 
lar mass, star-formation rate (SFR), dust extinction, and 
age. For MD5a, the clump without leaking LyC emission, 
we found: log(T) = 8.61 q g, log(Mass [Mq]) = 9.701q4q, 
log(SFR [Mo yr-i]) = 0.5ll°;ti E(B - V) = 0.02l°;J°, 
and log(Age [yr]) = 9.10lQ g 7 . In constrast, for MD5b, 
the clump with leaking LyC emission, we found: log(r) = 
8.21J-8, log(Mass [Mg]) = 8.691°-^^, log(SFR [M© yr'i]) 
= 1.3ltH?. -V) = 0.17t°;°3, and log(Age [yr]) = 

7.70tJ°'13. 

The confidence interval is bounded by log(Age [yr]) = 7.70 
(age = 50 Myr) because that is the minimum age we enforce upon 
the models. 
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The fits to MD5a and MD5b describe two very different 
stellar populations. MD5a has an old stellar population 
with significant stellar mass, a low star-formation rate, 
and very little reddening. MD5b, however, is young and 
low-mass, with a much higher star-formation rate and a 
larger E{B — V). While we have set a minim um LEG age 
limit of 50 Myr following iReddv et al.l (|2012ll , the best-fit 
model to MD5b without a minimum age requirement is 
even younger (10 Myr). Because of the differing stellar 
populations of MD5a and MD5b, it is possible that they 
are two distinct galaxies in the process of merging. It 
is also possible that they are simply two sub-regions of 
the same galaxy, one of which (MD5b) is undergoing a 
recent burst of star formation that has greatly increased 
its non-ionizing UV and LyC fluxes, making it more likely 
to detect LyC emission. 

In addition to modeling the stellar populations of 
MD5a and MD5b individually, we used FAST to model 
the best-fit stellar population to the combined photom¬ 
etry of both components. While such modeling does not 
represent a physically meaningful stellar populatiorF^. 
it facilitates a comparison to the results of ground-based 
LEG surveys, recreating the flux that would be measured 
by instruments lacking the high resolution of HST. For 
the combined photometry of MD5, we found: log(T) = 
8.2tJ;8, log(Mass [Mq]) = 9.52t0;|4, log(SFR [Mq yr'i]) 
= E{B -V) = and log(Age [yr]) = 

8.5 ±0.8. 

In order to compare the derived properties of MD5a, 
MD5b, and the combined photometry with those of typ¬ 
ical LBGs, we used FAST to fit SPS models to the re¬ 
mainder of the LBGs in our sample and examined their 
stellar masses, star-formation rates, dust extinction, and 
ages. We performed this analysis both for individual sub- 
arcsecond components of LBGs, for the combined fluxes 
from all components of each LBG (simulating ground- 
based studies that are free of contamination), and for a 
combined-flux sample that also includes foreground con¬ 
taminants in cases where the contaminant was indistin¬ 
guishable from the LBG in the ground-based LRIS imag¬ 
ing (a fair simulation of ground-based studies). Figure 
m shows histograms of these three LBG samples, along 
with the best-fit values for MD5a, MD5b, and the com¬ 
bined photometry for MD5. MD5a and MD5b clearly 
have distinct stellar populations from each other, and 
MD5b, the LyG-emitting component, stands out as be¬ 
ing among the youngest 10% galaxies in the BIST LBG 
sample. 

For all the LBGs in our sample, we now compare the 
stellar population fits to individual clumps with fits to 
the combined photometry. As expected, we find that 
model fits to individual galaxy components generally 
yield smaller stellar masses and SFRs than fits to the 
combined photometry. Also, individual sub-arcsecond 
components exhibit a wider range of reddening values 
than do galaxies with combined photometry, although 
the median reddening value is the same for both sam¬ 
ples. There is no significant difference in the median 
derived ages between the individual clumps and the com¬ 
bined photometry. Finally, we find that the addition of 

For example, the mass and SFR derived from the best-fit 
model to the combined photometry of MD5 are less than the 
summed individual masses and SFRs of MD5a and MD5b. 


foreground contaminants does not significantly alter the 
sample-averaged LBG properties, as foreground contami¬ 
nants close to LBGs are rare in the LBG sample without 
NB3420 detections. We conclude that occasional fore¬ 
ground contaminants superimposed upon LBGs do not 
greatly affect the stellar populations derived for galaxies 
in ground-based LBG surveys. 

In order to compare the properties of MD5 with those 
of a much larger parent sample of LBGs, we consider 
the set of 570 L BGs at 2.7 < z < 3 .7 from the ground- 
based survey of iReddv et al.l (|2012t) . In Figure [Ml we 
display parameters of the stellar population fit to MD5a, 
MD5b, and the combined p h otome try with respect to 
the LBGs from IReddv et al.l (120121) The stellar Popu¬ 
lation parameters we display for the IReddv et al.l (|2012D 
LBGs have been derived from the latest solar metallicity 
models of S. Chariot & G. Bruzual, using constant star 
formation histories and a min imum age lim i t of 5 0 Myr. 
The median parameters of the IReddv et al.l (I2012D LBGs 
are consistent with those of our combined-photometry 
LBG sample, although the percentage of young galaxies 
(<100 Myr) in our HST sample is less than that of the 
IReddv et al.l (I2012D sample 0 

To facilitate comparison with the IReddv et al.l (120121 1 
LBGs, we have re-fit the stellar populations of MD5a, 
MD5b, and the combine d photometry of M D5 using 
the methods described in IReddv et al.l (j2012tl , employ- 
ing CSF models, a minimum age of 50 Myr, and 
iCalzetti et al.l (j2000D extinction. These fit results for 
MD5 (plotted in Figure ITdl) are qualitatively similar to 
those from our original delayed-r fits. The data in Fig¬ 
ure [14] show that the fit to the combined photometry 
of MD5 is unremarkable when compared to the ground- 
based photometry of the IReddv et al.l (I2012D LBG sam¬ 
ple: it has a typical stellar mass and age, and slightly 
below-average values for SFR and E{B — V). It is only 
when MD5a and MD5b are fit separately that the young 
stellar population of MD5b becomes appar ent; an age of 
50 My r places it in the youngest third of the IR eddv et al.l 
(|2012l) LBG sample. Finally, we als o fit the phot o metry 
for MD5b using SMC extinction. IReddv et al.l (|201^ 
found that ~90% of LBGs with Calzetti-inferred ages of 
<100 Myr had older ages (>100 Myr) when fit using an 
SMC extinction curve. For MD5b, the best-fit model 
using SMC extinction had an age of 160 Myr, along 
with less reddening {E{B — V) = 0.05) than the best- 
fit Calzetti-attenuated model and qualitatively similar 
values for stellar mass and star-formation rate (log(Mass 
[Mf7)1) = 9.16, log(SFR [Mq yr“^]) = 0.95). IReddv et al.l 
(I2012D model their full LBG sam ple using a combina - 
tion of extinction curves, emoloving lCalzetti et al.l (I2000D 
attenuation for the majority of the sample, but using 
SMC extinction for galaxies with Calzetti-inferred ages 
younger than 100 Myr. Compared with the ages derived 
from these fits, the 160 Myr age estimated for MD5b is 
still in the youngest third of the sample. Table |3| sum¬ 
marizes the stellar population fits to MD5 with respect 
to those of the non-LyC-emitting galaxies in our HST 

The fact that we find few young galaxies (<100 Myr) among 
the i^5'T-imaged LBGs and LAEs in the HS1549 field may be due 
to statistical variation inherent to our small sample size (40 galaxies 
with ^7336^606^125-^^160 imaging), a peculiarity of the HS1549 field, 

S ssibly a property of all protoclusters (see, e.g., ISteidel et all 
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Figure 13. Histograms showing the stellar population parameters of MD5 calculated by FAST with respect to the distribution 
of parameters for all LBGs with HST Hsse V 606 di 25 Hi 60 imaging. MD5b (the LyC-emitting component) is indicated by the red 
vertical line, MD5a is indicated by the blue vertical line, and the combined photometry for MD5a and MD5b is indicated by 
the black dashed vertical line. Gray filled histograms indicate parameters derived from FAST fits to individual sub-arcsecond 
components of LBGs. Black hashed histograms indicate parameters derived from FAST fits to combined photometry, mimicking 
ground-based studies. Green open histograms also indicate parameters derived from fits to combined photometry, but include 
known foreground contaminants in cases where the contaminant was indistinguishable from the LBG in the ground-based LRIS 
imaging. The similarities between the green and black histograms show that the inclusion of foreground contaminants does not 
significantly alter the sample-averaged LBG properties. MD5b stands out as having an age in the youngest 10% of the LBG 
sample. 


Table 4 

Stellar Population Parameters 


Sample 

log(Stellar Mass) 

log(SFR) 

E{B - VY 

log(Age) 


[Mq] 

[Mq yr 1] 

[mag] 

[yr] 

MD5a^ 

9.70l°-j5 


0 n2+o-^o 

U.U^_Q Q2 

Q in+0.20 

y-J-U_ o.67 

MD5b (LyG detected)*’ 

8.69l°;« 

1 1 +0-09 

-‘-•'^-‘^-0.47 

0 1 7+0-03 
U.l<_o.o9 

7.70lo°® 

MD5 combined photometry^ 

Q C2+0-54 
y-O^_0.53 

1 1 q-|-0.58 
^•-*^'^-0.47 

0 10+0-^^ 
U.lU_o.io 

Q crj-|-0.80 

o.OU_o.8o 

LBG sub-arcsecond components'^ 

^•^’^-0.53 

0 75+0-76 
U./O_o.76 

0 1 0+0-20 
U.1U_Q 2Q 

q qr»-|"0.20 

o-of-0.50 

LBG combined photometry^ 

Q ■77+0.58 
‘ ' -0.26 

1 05+0-30 

i.uo_o.22 

0 10+0-^^ 
U.1U_Q QY 


LBG combined photometry -|- contaminants^ 

9-84tg:?i 

1 2Q+0-37 
4-^y-0.41 

0 10+0-^3 

U.1U_Q Q5 

q qr»-|-0.30 
O.OU_o.40 

2.7 < 2: < 3.7 LBGs from Reddv et al. (2012)^ 

9-78l°j^ 

1 65+0-40 

n 1 q+O-09 

b.lO_Q 10 

8.051°;®“ 


^ Derived from Ay, assuming k(V) = A:(5500) = 4.048 from ICalzetti et al.l 11200(11 1. 

^ Best fit values from FAST assuming a delayed-r star-formation history, a minimum age of 50 Myr, and 
ICalzetti et al.l II2000I 1 attenuation. Uncertainties quoted are 68% confidence intervals. 

Median values for LBG samples are quoted, along with values bracketing the inner 68% of the distributions. 
® Best-fit stellar population parameters for 2.7 < z < 3.7 L BGs from IR.ed dy et aP 11201211 . assuming a 
constant star-formation history, a minimum age of 50 Myr, and IGalzetti et al.l II2000I1 attenuation. Median 
values are quoted, along with values bracketing the inner 68% of the distributions. 
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LEG sample and those of the LBGs from iReddv et alJ 

dMl). 

As the model fit to MD5b has a young age and 
low stellar mass, two properties typical of LAEs (e.g., 
iGawiser et alJ . [2?)(T^ iGuaita et alJ . l2011fl . we also exam¬ 
ined MD5b with respect to the LAEs in our sample, none 
of which exhibited LyC detections. We modeled LAE 
stellar populations with FAST using the same methods 
as for LBGs, but in order to account for the young ages 
and lower metallicities associated with LAEs we set the 
minimum age to 10 Myr and fixed metallicity at 20% 
solar. Figure [15] shows the distribution of stellar mass, 
star-formation rate, dust extinction, and age for LAEs, 
along with values for MD5a, MD5b, and the combined 
photometry of MD5. MD5b has a stellar mass more typ¬ 
ical of the average LAE in our sample, but a higher than 
average SFR. Its age is still young compared to the LAE 
sample. 

While the young age of MD5b is shared by several ob¬ 
jects in our HST sample, none of these additional young 
objects exhibit LyG detections. In the LBG sample, 
two object subcomponents in addition to MD5b have 
ages less than 100 Myr, and five such components ex¬ 
ist in the LAE sample. If the young stellar population 
of MD5b is responsible for its LyC emission, then we 
might also expect LyC detections from other LBG and 
LAE components with similarly young ages. Both of the 
young objects in the LBG sample (M16a and MD34f) 
have small stellar masses (log(Mass [Mq]) ^ 9) and large 
SFRs (log(SFR [Mq yr“^]) > 1.6) like MD5b, but both 
of these objects are redder {E{B — E) ^ 0.4). If LyC 
emission is being produced copiously by the hot stars in 
these two galaxies, the additional dust extinction might 
be the reason we do not detect the LyC photons. We 
note that MD34f is actually associated with a NB3420 
detection (see Section [5.3.1|) . but its close proximity to 
a foreground contaminant makes it impossible to dis¬ 
tinguish between emission from MD34f and the fore¬ 
ground contaminant in the seeing-limited NB3420 im¬ 
age. The LAE sample presents several additional young 
galaxy components that are undetected in the NB3420 
image. These objects have reddening values similar to 
MD5b {E{B — V) ^ 0.2), lower than those of the young 
LBGs. One possible explanation for the lack of LyC de¬ 
tections in the LAE sample is that MD5b is located along 
a fortuitously clear IGM sightline, and these LAEs are 
not. Another possibility is that the LAEs are simply 
too faint to be detected in our LyC imaging. The Eeoe 
magnitude of MD5b is meoe = 25.9, and its LyC detec¬ 
tion is near the edge of our detection limit. The Eeoe 
magnitudes of the young LAEs in our sample (ages < 
100 Myr) are much fainter on average, and range from 
26.29 < rriLyC < 29.14, with a median of 27.69. If these 
objects have the same ratio of ionizing to non-ionizing 
radiation as MD5b, the LyC magnitudes of these objects 
would range from 27.80 < meoe < 30.64, with a median 
of 29.20. Such magnitudes are well below the detection 
limit of the NB3420 filter used for LyC imaging (27.3 
mag), and thus these objects would not have been de¬ 
tected in the LyC. 

6.4. No Evidence for AGN 

IVanzella et al.l (|20I5ll discuss the possibility that some 
portion of the leaking ionizing radiation from the LyC- 


emitters lonl and Ion2 are due to contributions from 
low-luminosity AGNs within these galaxies. We consid¬ 
ered this possibility in MD5b, but find no evidence for a 
low-luminosity AGN. There is no significant variability 
detected between the E—band magnitude of MD5 mea¬ 
sured by LRIS in 2007 (my = 24.96 ± O.II) and the 
combined HST Eeoe magnitudes of MD5a and MD5b 
measured in 2013 (meoe = 25.11 ±0.08). Additionally, 
the available spectra of MD5 do not show any high ion¬ 
ization emission line features (e.g., NV A1240 emission), 
although our spectrum does not cover CIV A1550 or Hell 
A1640, nor are we able to examine the OIII / H/3 ratio 
because H/3 falls on a sky line. We therefore find no ev¬ 
idence for AGN activity in MD5b with the information 
available. 

7. DISCUSSION 

In this section, we consider the broader implications 
of our single robust detection of LyC emission in the 
HS1549 field. We discuss the actual rate of foreground 
contamination with respect to predictions from the con¬ 
tamination simulations of M13, and the LyC emission 
properties of MD5b, including its intrinsic ionizing to 
non-ionizing UV flux-density ratio and the implied LyC 
escape fraction. With the removal of all foreground con¬ 
taminants from the M13 sample, we then obtain a revised 
estimate for the ionizing emissivity due to star-forming 
galaxies at z = 2.85. Finally, we discuss the prospects for 
future direct searches for LyC radiation in high-redshift 
galaxies. 

7.1. Comparison to Previous Contamination Estimates 

We wish to address the question of whether or not the 
simulations of foreground contamination from MI3 ac¬ 
curately predicted the number of contaminants in the 
sample. Th e simu lations (described in detail in MI3; 
iNestor et al.l . I20I3I1 employed the surface density of ob¬ 
jects in the NB3420 image and the offset of each NB3420 
detection to estimate the number of foreground contami¬ 
nants and a contamination-corrected NB3420 magnitude. 
The simulations predicted that 1.5 ± 1.0 out of 4 LBGs 
and 4.3 ± 1.3 out of 7 LAEs in the main sample were 
real LyC-emitters. The prediction for the LBG sam¬ 
ple has held out, as one LBG (MD5) remains a strong 
candidate for LyC emission. The prediction for the 
LAE sample, however, was too high. We were able to 
obtain t/aaeVeoe >^ 125^^160 imaging for 6 out of 7 LAEs 
with LyC detections, yet none have proven to be true 
sources of LyC emission. We note that 2 of the 7 LAEs 
had misidentified redshifts {lae2436 and lae7180), and 
should not have made it into this sample in the first 
place. Given the small number of galaxies with true 
LyC detections, contaminated NB3420 detections from 
even one or two objects with misidentified redshifts may 
introduce a non-negligible bias that is not taken into ac¬ 
count in the contamination simulations of M13. Accord¬ 
ingly, we re-ran the same contamination simulations, but 
only considering the 4 spectroscopically-confirmed LAEs 
from the main sample of M13 for which we acquired 
U 33 eVGoeJi 25 Hieo imaging. The revised simulations pre¬ 
dict 2.1 ± 1.0 out of 4 galaxies to be real LyC-emitters. 
As none of the four NB3420 detections proved to be real, 
the prediction is still too high by ~ 2cr. With the small 
sample size of 4 objects, however, such variations may 
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Figure 14. Histograms showing the stellar populat ion paramete r s of M D5 calculated by FAST with respect to the distribution 
of parameters for the 2.7 < 2 < 3.7 LEG sample of IReddv et all (I2012li . MD5b (the LyC-emitting component) is indicated by 
the red vertical line, MD5a is indicated by the blue vertical line, and the combined photometry for MD5a and MD5b is indicated 
by the black dashe d vertical line. All parameters plotted are derived using the stellar population fitting methods descr i bed in 
IReddv et al.l (120121 1. with assumptions of constant star formation histories, a minimum age of 50 Myr, and ICalzetti et ^ (1200011 
dust attenuation. While the propertie s of th e composite object MD5 (containing both MD5a and MD5b) do not stand out 
among other LBGs in the IReddv et al.l (1201211 sample, the LyC-emitting component MD5b stands out as being in the youngest 
third of the sample. 


be expected. We conclude that the contamination sim¬ 
ulations still serve as a useful, though blunt, tool for 
evaluating the likelihood of foreground contamination. 

7.2. Ionizing to Non-ionizing UV Flux-density Ratios 

One of the intriguing findings reported by all 
three ground - based LyC studies o f gala xy protoclusters 
(jlwata et al.L 120091 : INestor et al.l . 120111 : M13) was the 
high apparent ratio of ionizing to non-ionizing radia¬ 
tion in many of the candidate LyC-emitters, which ap¬ 
peared to be in conflict with results from standard SPS 
models. Seyer al models haye bee n proposed to address 
this question. Ilwata et al.l (|2009ll proposed a top-heavy 
IMF. iNestor et al.l (|2013ll investigated the intrinsic non¬ 
ionizing to ionizing UV flux-density ratios for two sets 
of stellar population models at varyi ng ages and metal- 
licites . These authors examined both iBruzual fc Charlotl 
(j2003D models and BPASS models, which include a more 
detailed treatment of stellar binaries and Wolf-Rayet 
stars and result in bluer galaxy spectra. In this work, 
we also consider S B99 models with impr oved treatment 
of stellar rotation (|Leitherer et al.L 120141) , as rapidly ro¬ 
tating stars exhibit bluer spectra as well. 

However, our new observations, along with all pre¬ 
vious follow-up work aimed at investi gating contami- 
nation among candidate LyC-emitters (jVanzella et al.l . 


120121:ISiana et al.l . 120151) . have ruled out all objects with 
high apparent ratios of ionizing to non-ionizing radia¬ 
tion as contaminants. In our sample of 16 galaxies with 
putative LyC emission, 10 had anomalous flux-density 
ratios of {Fuv/FLyc)obs < 2.0. However, the results 
of our analysis show that MD5 remains the only robust 
candidate for LyC emission. The flux-density ratio of 
MD5, as calculated from the ground-based NB3420 and 
U—band imaging of M13, is Fuv/FLyC = 5.9 ± 2.0. Us¬ 
ing our HST U 336 and Vqob imagin eFI to revise this 
flux-density ratio to only include MD5b, the compo¬ 
nent associated with the ionizing radiation, we obtain 
a value of Fuv/FLyC = 4.0 ± 2.0. At the redshift 
of MD5 (z = 3.14), the maximum IGM transmission 
through the U 33 Q filter is ^ 60 — 7 0% (depending on 
the absorbe r statis tics used; see, e.g.. lRudie et al.l . 120131 : 
llnoue et al.l . 120141) . although IGM transmission varies 
greatly along the line of sight. These maximum transmis¬ 
sion values give rise to IGM-corrected flux-density ratios 
of Fjjv/FLyC = 2.4 — 2.8, demonstrating that the flux- 
density ratio of MD5b is consistent with expectations 
from BPASS models between the ages of ~ 10 — 50 Myr 

We use 1/336 rather than NB3420 to represent flux in the LyC 
because the matched-PSF, matched-aperture photometry between 
the f/336 and Veoe images allows the most accurate calculation of 
Fuv/FhyC- 
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Figure 15. Histograms showing the stellar population parameters of MD5 calculated by FAST with respect to the distribution 
of parameters for all LAEs with HST UzmV&o&J\ 25 H\Qo imaging. MD5b (the LyC-emitting component) is indicated by the red 
vertical line, MD5a is indicated by the blue vertical line, and the combined photometry for MD5a and MD5b is indicated by 
the black dashed vertical line. Histogram colors are as in Figure [T^l The bottom-right plot shows several LAEs with log(age) 
< 8, similar to MD5, raising the question of whether or not these objects — which do not have NB3420 detections — may also 
emit LyC radiation. However, as discussed in Section [6.31 these LAEs are much fainter than MD5. If their observed ionizing 
to non-ionizing UV flux density ratios were equivalent to that of MD5, then they would be too faint to detect in the NB3420 
image. 


(see Table 7 of INestor et al.l . [201311 in the lack of signifi¬ 
cant IGM absorption. 

In similar work inv estigating c o ntami nation in 
candidate LyC-emitters, iSiana etlill (|2015D obtained 
spatially-resolved spectroscopy of the five LEG candi¬ 
dates for LyC emission presented in iNestor et al.l (120111) . 
While unable to confirm any candidates, their data 
showed that both galaxies in their sample with anoma¬ 
lously high apparent ratios of ionizing to non-ionizing 
radiation (MD32, aug96ml6) are contaminated by lower- 
redshift objects. Add itionally, neither of the L yC- 
emitters lonl or Ion2 (|Vanzella et al.l . I2012L1201511 ex¬ 
hibit anomalously high ratios of ionizing to non-ionizing 
radiation. While the two dozen LyC-emitter candidates 
in the literature with extensive follow-up data do not 
comprise the entire sample of high-redshift candidate 
LyC-emitters, the fact that all candidates with extreme 
{Fuv/FLyc)obs ratios have proven to be contaminants 
may indicate that the anomalously high ratios of ionizing 
to non-ionizing radiation originally inferred are simply a 
result of foreground contamination. 

7.3. Escape Fraction for MD5 

We can estimate the relative and absolute escape 
fractions of MD5b using the intrinsic luminosity- 
density ratio {Luv/LLyc)intr, observed flux-density ra¬ 


tio {Fjjv / FLyc)obs, and IGM transmission factor (tjoM) 
with the following equations: 

nLyC ^ {Luv/LLyc)^ntr (-in 

esc.rel {FuV/FLyc)obs ^ ^IGM 

pLyC _ pLyC fUV in\ 

J esc.abs Jesc.rel J esc 5 V^/ 

The intrinsic luminosity-density ratio is highly un¬ 
certain, and varies significantly with the stellar popu¬ 
lation synthesis models used, as well as with the age 
and metallicity of t he stellar population (see Table 7 of 
iNestor et al.l . I20I3D . Here, we bracket the full range of 
possible values for the LyC escape fraction of MD5b by 
using the lowest value of {Ljjv / LLyc)intr = 2.1 (from 
BPASS models of age lOMyr with 20% solar metallicity) 
and a maximum IGM transmission of 70%. From these 
values, we obtain fesc,rei = 75%, which must be inter¬ 
preted as a lower limit such that the range of allowed rel¬ 
ative escape fractions is — 75 — 100%. We use the 

value = 0.19 to estimate the UV escape fraction at 
1500A, which is calculated from the E(B — V)oi the best- 
ht FA ST model to MD5b, assuming the iCalzetti et al.l 
(j200flH attenuation curve. We thus obtain an absolute 
escape fraction of f^^abs = 14 — 19%. As the flux 
has photometric errors of roughly fifty percent, the un- 
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Table 5 

Contributions to the Ionizing Background. 


FuvlFLyC° Magnitude range'" <^Lyc'^ 


LyC Detections: MD5 only® 


(i) 

LBG 

140+'^ 

Mab < -19-7 

0.8 ±0.3 

(ii) 

LAE 

> 14 

-19.7 < Mab < -17.7 

< 1.7 

(iii) 

LBG 

> 14 

-19.7 < Mab < -17.7 

< 7.3 

(iv) 

LBG 

140+™ 

Mab < -17.7 

1.5 ±0.6 

(v) 

LAE 

> 14 

Mab < -17.7 

< 3.6 


Total (lum.-dep.)^ 
Total (LAE-dep.)^ 


Mab < -17.7 
Mab < -17.7 

0.8 ± 3.7 
1.2 ± 1.9 


LyC Detections: 

MD5, D24, lae46S0^ 


(i) 

LBG 

'^-16 

Mab < -19.7 

1.5 ±0.4 

(ii) 

LAE 

71 

'-*--18 

-19.7 < Mab < -17.7 

0.3±0.1 

(iii) 

LBG 

71 

'-*--18 

-19.7 < Mab < -17.7 

1.4 ±0.5 

(iv) 

LBG 

74+24 

'^-16 

Mab < -17.7 

2.9 ±0.8 

(v) 

LAE 

-jr-j +34 
'-*--18 

Mab < -17.7 

0.7±0.2 


Total (lum.-dep.)^ 
Total (LAE-dep.)s 


Mab < -17.7 
Mab < -17.7 

3.0 ± 0.9 
2.9 ± 0.8 


^ Luminosity function parameters are identical to those in M13. 

^ Sample average flux-density ratio corrected for IGM absorption, described in Section 17.41 
® Magnitude range over which the first moment of the luminosity function is determined. 

= —19.7 and —17.7 correspond to 0.34L* and 0.06L*, respectively. 

^ Comoving specific emissivity of ionizing radiation in units of 10^^ ergs Hz“^ Mpc“^. 

® As we do not have the full suite of f/aaeVgoe•^125-^^160 imaging for D24 and lae4680, we cannot 
determine whether their NB3420 detections are true LyC emission or foreground contamina¬ 
tion. Therefore, we perform two calculations of the emissivity in order to determine the full 
range of its possible values. In the upper portion of this table we assume that both D24 and 
lae4680 are foreground contaminants, and that MD5 is the only true LyC detection. In the 
bottom portion of the table, we assume that MD5, D24, and lae4680 are all true LyC-emitters. 
^ Total for the luminosity-dependent model, determined by su mmi ng rows (i) and (iii). Limits 
are taken into account using the method described in Section \7a\ 

® Total for the LAE-dependent model, determined by summing 0.7 7x r ow (iv) and row (v). 
Limits are taken into account using the method described in Section 17.41 


certainty in f^sc is also at minimum fifty percent. 

7.4. Revised LyC Emissivity for Star-forming Galaxies 
at z = 2.85 

Here, we present a revised calculation of the emissivity 
of ionizing photons at z = 2.85 based on the analysis 
of the HST data in the HS1549 field. We estimate the 
comoving specific emissivity as 

eLyC=(§^) ' (3) 

follow ing the assumptions of M13 and iNestor et al.l 
(|2013f) . In this expression, L is the non-ionizing UV 
luminosity, $ is the non-ionizing UV luminosity func¬ 
tion, and [Fuv/FLyc)corr is the average flux-density ra¬ 
tio of non-ionizing to ionizing UV radiation for the en¬ 
tire galaxy sample, corrected for the mean IGM atten¬ 
uation in the LyC spectral regiorP^. We perform this 
emissivity calculation separately for the main sample of 
spectroscopically confirmed LBGs and LA Es from M13 
using the UV lumi nosity functions fr om iR.eddv et al.l 
200^ for LBGs and lOuchi et al.l (1200811 for LAEs), and 
combine the LEG and LAE emissivities to obtain a to¬ 
tal emissivity for star-forming galaxies. As in M13, we 
use the LRIS U—band to represent non-ionizing UV flux 
and NB3420 to represent LyC flux. The difference be- 

To correct for absorption of LyC photons by neutral hydro¬ 
gen in the IGM, we use the sample-averaged transmission values 
calculated in M13: = 0.44 ± 0.03 and = 0.35 ± 0.04. 


tween our calculation and that of M13 lies in our estima¬ 
tion of the average flux-density ratio. Rather than esti¬ 
mating the average amount of foreground contamination 
from simulations, we instead know exactly which galax¬ 
ies are contaminated based on the HST data. There 
were only two NB3420-detected galaxies in the M13 spec¬ 
troscopic sample for which we were unable to acquire 
f^ 336 b 606 >^i 25 -fli 60 imaging (D24 and lae4680), and for 
these objects we could not evaluate whether or not the 
NB3420 detections are due to foreground contamination. 
We thus calculate the emissivity twice in order to quote 
the full range of possible values: in one calculation we 
assume that MD5 is the only true LyC detection, and 
in the other calculation we assume that MD5, D24, and 
lae4680 are all true LyC-emitters. In addition to us¬ 
ing the HST data to remove the NB3420 flux of fore¬ 
ground contaminants, we also use these measurements 
to estimate the percentage of contaminated flux in the 
non-ionizing UV. All objects with foreground contami¬ 
nants identified through the HST imaging are blended 
in the LRIS V imaging, and it is impossible to isolate 
the uncontaminated z ~ 2.85 flux in the LRIS V im¬ 
age. Thus, for each contaminated object, we decrease its 
LRIS U—band flux to match the fraction of uncontami¬ 
nated Veoe flux in the HST imaging. For objects that do 
not have HST C/ 336 V 606 >^i 25 ^Ii 60 imaging and are unde¬ 
tected in NB3420, we decrease their LRIS U—band flux 
to match the average fraction of uncontaminated Veoe 
flux in the full sample of HST-iva&ged galaxies with¬ 
out NB3420 detections (99% for LBGs, 91% for LAEs). 
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Finally, we use the same sample-averaged IGM correc¬ 
tion to compute {Fuv/FLyc)corr as descri bed in M13, 
emplo ying statistics of HI absorbers from iRudie et alJ 
(j20Rlfl . We note that the clustering of Lyman limit sys¬ 
tems is not taken into account in these absorber statis¬ 
tics, and thus the true mean IGM transmission may be 
slig htly higher than the va lues presented in MIS (see, 
e.g., iProchaska et alJ . I2Q14I1 . 

In order to estimate the total contribution of star¬ 
forming galaxies to the ionizing emissivity at z ^ 2.85, 
we estimate the emissivity due to LBGs and LAEs sep¬ 
arately a£d__use_two different models (described in de¬ 
tail in iNestor et al.l . 1201311 to combine these values into 
a total emissivity of star-forming galaxies. In the first 
model, which we refer to as the luminosity-dependent 
model, the Ffjv/FLyC values for LAEs are assumed to 
represent those for star-forming galaxies with faint UV 
continuum magnitudes (0.06L* < L < 0.34L*, corre¬ 
sponding to 25.5 <V< 27.5) while the Fuv/FLyc val¬ 
ues for LBGs represent those for brighter star-forming 
galaxies {L > 0.34L*). The second model, referred to as 
the LAE-dependent model, considers the case that LAEs 
are not simply faint LBG-analogs, but that LBGs and 
LAEs are actually distinct populations of star-forming 
galaxies with systematically different Fjjv/FLyC values 
on average. In this model, LAEs are assumed to com - 
prise 23% of the LBG population (INestor et al.l . 1201311 . 
galaxies identified both as LBGs and LAEs are treated 
as LAEs, and the luminosity function is integrated over 
the full luminosity range (0.06L* < L < oo) for both 
LBGs and LAEs. 

In Table [SJ we summarize the contributions to CLyC as 
determined from galaxies in the HS1549 field. By consid¬ 
ering MD5 as the only real LyG detection in the entire 
galaxy sample, we obtain values of the average UV flux- 
density ratio for LBGs and LAEs the HS1549 field to be 
{FuvjFLyct^,^ = and {FuvIFLyc)tif > W- 

If we include the NB3420 detections for D24 and lae4680 
as well, these values become {Fuv/FLyc)cm-r = 
and {Fuv/FLyc)c<^r = The uncertainties in the 

flux-density ratios are dominated by the NB3420 photo¬ 
metric errors, and are large because of significant uncer¬ 
tainty in the average NB3420 flux due to our very few 
NB3420 detections. We estimated the uncertainties in 
the flux-density ratios using a Monte Garlo simulation. 
The simulation calculates a distribution of average val¬ 
ues for {Fuv!Fuyc)corr based on random realizations of 
our photometric data and IGM correction within their 
Gaussian uncertainties. As the resulting distribution 
is positively skewed, we quote error bars that bracket 
the inner 68% of the distribution. Based on the val¬ 
ues for {Fuv/FLyc)corr, we infer revised values of the 
comoving specific emissivity: considering NB3420 emis¬ 
sion from MD5 only (or from MD5, D24, and lae4680) 
we obtain euyC = 0.8 ± 3.7 (3.0 ± 0.9) x 10^'^ ergs s“^ 
Hz“^ Mpc“^ for the luminosity-dependent model and 
euyC = 1-2 ± 1.9 (2.9 ± 0.8) x lO^^ ergs s"! Hz"! Mpc-^ 
for the LAE-dependent model. The uncertainties in euyC 
reflect only uncertainties in {Fuv/FLyc)corr, which dom¬ 
inate over uncertainties in the luminosity function. Each 
total emissivity value represents the linear combination 
of the LBG and LAE contributions to the emissivity, 
which are characterized by distinct {Fuv/FLyc)corr val¬ 


ues and integrals over the non-ionizing UV luminosity 
functions. To determine the error bar for each contri¬ 
bution to the emissivity, we divided the relevant non¬ 
ionizing luminosity function integral by the correspond¬ 
ing random distribution of {Fuv/FLyc)corr values de¬ 
scribed above. We then randomly drew linear combi¬ 
nations of emissivities (a simple sum for the luminosity 
dependent model, and a weighted linear combination for 
the LAE-dependent model; see Table 5, notes f and g) 
from the relevant distributions. To combine the LBG 
and LAE emissivities in cases where the value of the 
LAE emissivity is only an upper limit, we treated the 
LAE emissivity as a normally distributed random vari¬ 
able centered on zero, with a standard deviation equal 
to its Icr limit. The resulting distribution of emissivities 
was not skewed, so we report the standard deviation of 
the distribution as the uncertainty in the emissivity. 

The revised values of euyC are much lower than those 
computed in MI3: euyC = 15-0 ± 6.7 x 10^^ ergs s“^ 
Hz“^ Mpc“^ for the luminosity-dependent model and 
^LyC = 8.8 ± 3.5 X 10^^ ergs s“^ Hz“^ Mpc“^ for the 
LAE-dependent model. The lower emissivity value cal¬ 
culated in the current work is much more compatible 
with the total ionizing emissivity at 2 ; = 2.85, estimated 
in MI3 from measurements of the Lyg-forest opacity 
(jBolton fc HaehneI9 . l2007l:lFaucher-Giguere et aLll^OO^ 
to be e*'£yu; ^ 5 — 10 x 10^"^ ergs s“^ Hz“^ Mpc“^. 
As measurements of the contribution of QSOs to the 
ionizing background at z = 2.85 range from ~ 

1.5 X 10^"* ergs s“^ Hz“^ Mpc (|Cowie et al.L 120(391) to 
^Lu (7 5-5 X 10^"^ ergs s“^ Hz“^ Mpc (|HoDkins et al.l . 

[2001 . our data indicate (with large uncertainties) that 
star-forming galaxies provide roughly the same contribu¬ 
tion as QSOs to the ionizing background at this redshift. 


7.5. The Future of LyC Surveys 

The results from this work suggest that identifying true 
LyC-emitters at high redshift requires an extremely large 
parent sample of galaxies and/or significantly deeper 
LyC observations. With only one confirmed detection 
out of 49 LBGs, and zero confirmed detections among 
the 91 LAEs, the detection rate of LyC emitters at high 
redshift is very small. While several interesting meth¬ 
ods of indirectly identifying LyC-emitting galaxies have 
been proposed — such as assessing the shape of the 
Lya emission line (jVerhamme et all 1201 41). observing re¬ 
duced flux in nebular emission lines ([Zackrisson et aB . 
I 2 OI 3 D . and observing residual flux in the cores of satu¬ 
rated low-ionization a bsorption lines (iJones et al.l . 120131 : 
iHeckman et al.l . 120111 ) — it will not be possible to verify 
the validity these indirect methods without first obtain¬ 
ing a sample of galaxies with robust detections of LyC 
emission. 

Given their low detection rate in the HS1549 field, it 
will likely be very difficult to amass a statistical sample 
of LyC emitters at 2 < z < 4 without a dedicated survey. 
As the process of identifying and verifying LyC emission 
in this redshift range involves several steps, we outline 
here what we consider to be the most efficient method 
for doing so. 

First, there is the question of efficient target¬ 
ing. Because the average surface density of LBGs 
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down to i? = 2 5.5 is rough l y 1 — 2 ga laxies per 
square arcminute (iSteidel et al.l . I1999L l2004f) . the pro¬ 
cess of observing LBGs for LyC emission is greatly 
streamlined by observing galaxy protoclusters, which 
have an increased density of objects at a particu¬ 
lar redshift. Several galaxy protoclusters have al¬ 
ready bee n identified in t h e lite r ature at 2 < z < 4 
(see, e.g.. IVenemans et all 120071: [Kodama et all 12007 : 


Hatch et ah . 20HT Cucciati et all 2014 Lemaux et al 


2014 IShimal^wa et al.l . 12014 IDiener et al.l . I2015|b At 

the rate of one LyC detection per protocluster, observa¬ 
tions of at least ~ 10 protoclusters would be necessary to 
obtain a sample large enough to investigate systematic 
differences between LyC-leakers and non-leakers. While 
the environment in the ICM surrounding protoclusters 
may not be typical of the universe as a whole, it is un¬ 
likely to affect the escape of ionizing photons through 
the ISM of LyC-emitting galaxies. The factor that is less 
well-constrained in protocluster environments is the es¬ 
timate of the sample-averaged ICM transmission, which 
enters into the calculation of the global ionizing emis- 
sivity. However, galaxies in protoclusters can still be 
very useful for studying the multiwavelength properties 
of LyC emitters. An additional potential problem with 
targeting protoclusters for LyC studies is that protoclus¬ 
ters may be composed of galaxies with older stellar pop- 
ulations on average (as we found in this work; also, see 
ISteidel et al.l . I2005f) . If LyC emission is primarily emit¬ 
ted from galaxies with younger stellar populations (as 
suggested by the LyC detection for MD5b), then LyC 
emitters may be less common in protoclusters. 

Next, there is the question of the required observa¬ 
tions. Spectroscopic redshifts must be measured for a 
large sample of protocluster galaxies, as photometric red- 
shifts are not sufficiently precise to determine whether or 
not apparent LyC emission originates below 912A. Mul¬ 
tiwavelength HST imaging should then be obtained to 
probe the SED shapes and LyC emission for all galaxy 
components near the high-redshift targets (as in this 
work). To allow for direct LyC imaging, the protocluster 
identified must be at a redshift where one of the cur¬ 
rently available HST filters probes the LyC spectral re¬ 
gion just bluewards of the Lyman limit (such as F336IE 
for z > 3.06, or F275W for z > 2.38). To obtain useful 
limits on the amount of escaping ionizing radiation from 
the faintest galaxies, LyC magnitudes must be probed 
several times fainter than their non-ionizing UV magni¬ 
tudes. The faintest galaxies in our sample, LAEs, have 
Veoe magnitudes ranging from 26.08 < m-eoe < 29.14 
with a median of 27.62. In order to measure ionizing 
to non-ionizing flux-density ratios equivalent to that ob¬ 
served for MD5b {Fjjv / PhyC 4) for the faintest LAE 
(iTieoe = 29.14), the required LyC observations must 
reach a depth of ~ 30.7 magnitudes. Assuming object 
sizes close to the PSE size, this requires imaging roughly 
4 times more sensitive than the t /336 observations in 
this current work. Larger objects, comparable in size 
to MD5b, would require imaging ~ 20 times deeper. In 
practice, the best way to measure the average Fjjv / ^LyC 
ratio for the faintest galaxies may be with stacked LyC 
observations, or with t he next-generat i on UV space tele¬ 
scope (e.g., ATLAST: iPostman et al.l . 1200911 . Finally, if 
the morphology of the candidates for LyC emission are 


complex, the last step would be to obtain high-resolution 
spectroscopic follow-up of the LyC-emitting component 
of the galaxy. 

This plan is streamlined relative to the process we 
have followed thus far because it skips the time-intensive 
ground-based LyC imaging and analysis. The main ben¬ 
efit of the ground-based LyC imaging was that we were 
able to design a custom, narrowband filter for the exact 
redshift of the protocluster. Narrowband filters placed 
just bluewards of the Lyman limit are the least affected 
by ICM absorption, and probe LyC emission at wave¬ 
lengths where the LyC photons are most likely to ion¬ 
ize hydrogen. Even so, these benefits do not outweigh 
the cost in time and resources if HST filters are already 
available at the correct wavelengths to probe LyC emis¬ 
sion. For cluster redshifts where HST filters are not 
available for LyC imaging, but the Lyman limit falls 
above the atmospheric cut-off, it would be possible to 
obtain ground-based narrowband LyC imaging first at 
z ^ 3 with, e.g., Keck/LRIS. The roughly two dozen 
LyC-emitter candidates that would fall within a single 
Keck/LRIS pointing (most of which would be contam¬ 
inants) could be followed up individually with an AO- 
assisted integral fiel d spectrograph such as Keck/OSIRIS 
(|Larkin et al.l . [200^ . In all future LyC searches, it is im¬ 
perative to obtain high-resolution imaging and redshift 
confirmation of each galaxy sub-arcsecond component as¬ 
sociated with apparent LyC emission in order to rule out 
foreground contamination. 

8 . SUMMARY 

In MI3, we identified 30 candidates for LyC emission 
via detection in the Keck/LRIS NB3420 filter: 5 LBCs 
and 7 LAEs spectroscopically confirmed at z > 2.82 (the 
main sample), 10 photometric LAE candidates, and 8 
spectroscopically-confirmed LAEs not part of the main 
sample. In this current work, we have presented follow¬ 
up HST t/ 336 V 606 ^i 25 l?i 60 observations of 16 of these 
objects: 4/5 LBCs and 6/7 LAEs in the main sample, 
4/10 photometric LAE candidates, and 2/8 LAEs outside 
of the main sample. 

In our high-resolution HST imaging, all of the candi¬ 
dates for LyC emission exhibit significant substructure. 
We have thus used the HST imaging to obtain photomet¬ 
ric redshifts of each galaxy sub-arcsecond component in 
order to determine if the source of the NB3420 emission is 
truly at z > 2.82, or if it is from a lower-redshift contam¬ 
inant. Of the sixteen galaxies with NB3420 detections 
imaged in t/ 336 V 606 </i 25 //i 60 : nine were located near fore¬ 
ground contaminants responsible for the NB3420 emis¬ 
sion. Two objects had incorrect redshifts assigned to 
them, and thus the NB3420 emission was also from a 
low-redshift galaxy. Four objects, all LAEs without spec¬ 
troscopic confirmation, exhibited ambiguous SED shapes 
consistent with both z ^ 2.85 galaxies and foreground 
contaminants. Lack of spectroscopic redshifts for these 
objects, combined with their ambiguous SED shapes, 
makes it impossible to verify their NB3420 detections as 
true LyC emission. In the end, only one robust candidate 
for LyC emission remained: the LBC, MD5. 

MD5 has a spectroscopic redshift of z = 3.14, mea¬ 
sured from LRIS and MOSFIRE spectra containing Lya 
and [OHI] emission lines, along with multiple coincident 
interstellar absorption features, blueshifted with respect 
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to the systematic redshift and indicative of an outflow 
{Avis ^ 280 km s“^). These spectra show no spurious 
emission or absorption features indicating a foreground 
contaminant. In the HST imaging, it becomes appar¬ 
ent that MD5 is composed of two clumps (MD5a and 
MD5b), which may either be two components of the 
same galaxy or two separate galaxies in the process of 
merging. MD5b is associated with the NB3420 detec¬ 
tion. The best-fit stellar population synthesis model to 
MD5b indicates that while values for its stellar mass and 
reddening are typical of LBGs, it has a young stellar 
population (<50 Myr) and a high SFR (20 Mq yr“^) for 
such a low-mass object (5 x 10® Mq). This age places 
MD5b in the youngest 10% of the HST sample, and 
in the youngest third of typical LBGs. While MD5b 
is young compared to the full sample regardless of the 
fitting methods, we caution that the exact value for the 
best-fit age depends significantly on the dust attenuation 
curve and star-formation history assumed. 

The observed non-ionizing to ionizing UV flux-density 
ratio of {Ftjv/FLyc)obs = 4.0 ± 2.0 of MD5b is consis¬ 
tent with predictions of the intrinsic flux-density ratio for 
galaxies of 10—50 Myr from stellar population synthesis 
models (|Nestor et all 1201^ . in the absence of significant 
IGM absorption. With the assumption of maximum IGM 
transmission at the redshift of MD5b (t/cM = 70%), the 
observed flux-density ratio results in a relative escape 
fraction of = 75 — 100% and an absolute escape 

fraction of f^}^abs = 14 — 19%. We also note that the 
emission in the t /336 filter at the location of MD5b, which 
probes LyG emission at z = 3.14, shows no spatial off¬ 
set from the Veoe emission, supporting the interpretation 
that the ionizing photons are escaping either isotropi¬ 
cally, or through a hole in the ISM directly along our 
line of sight. 

MD5b stands out as one of the youngest galaxies in 
our HST-imaged sample. The best-fit model to MD5b 
had the youngest age available (50 Myr), and the best- 
fit becomes even younger (10 Myr) without the mini¬ 
mum age requirement. While there are an additional 
two LBGs and five LAEs in our sample with ages <100 
Myr, none of these galaxies exhibit NB3420 detections. 
The two young LBGs exhibit increased dust extinction 
{E{B — V)^ 0.4), which may be the reason we do not 
detect the LyG photons. The young LAEs have less dust 
extinction {E{B — V) ^ 0.2), but have very faint UV con¬ 
tinuum magnitudes (median meoe = 27.7). If the young 
LAEs have the same ratio of non-ionizing to ionizing radi¬ 
ation as MD5b, then their median LyC magnitude would 
be rriLyC = 29.20, well below the detection limit of the 
NB3420 filter used for LyG imaging (27.3 mag). 

After eliminating foreground contaminants from our 
sample, we obtained a revised estimate for the comoving 
specific ionizing emissivity {eLyc) at z = 2.85. We cal¬ 
culated the emissivity associated with LBGs and LAEs 
separately, and combined these va lues using two differ¬ 
ent models, described in Section 17.41 If we consider 
MD5 as the only galaxy with a LyG detection, we obtain 
^LyC = 0.8 ± 3.7 X 10^^ ergs s“^ Hz“^ Mpc“® for the 
luminosity-dependent model and esyC = 1-2 ± 1.9 x 10^^ 
ergs s“^ Hz“^ Mpc“® for the LAE-dependent model. 
If we also add in as true LyG detections the only two 
galaxies (D24 and lae4680) for which we were unable 


to obtain sufficient HST data to evaluate contamina¬ 
tion, we obtain esyC = 3.0 ± 0.9 x 10^'^ ergs s“^ Hz“^ 
Mpc“® for the luminosity-dependent model and esyC = 
2.9 ± 0.8 X 10^^ ergs s“^ Hz“^ Mpc“® for the LAE- 
dependent model. These revised values of esyC are much 
lower than those computed in M13, and much more com¬ 
patible with the total ionizing emissivity at z = 2.85 
{esyc = 5.6 ± 1.6 X lO^^ ergs s"! Hz"! Mpc"®; M13). 
Within the large photometric uncertainties, and uncer¬ 
tainties due to the small dynamic range in which we can 
probe LyC emission, our data indicate that star-forming 
galaxies provide roughly the same contribution as QSOs 
to the ionizing background at this redshift. 

Overall, the rate of foreground contamination for ap¬ 
parent LyC leakers in our z ^ 2.85 sample was very high. 
While the single detection in the LBG sample is con¬ 
sistent with contamination expectations from M13, the 
contamination rate in the LAE sample was higher than 
predicted by our simulations. With this work we have 
shown that ground-based LyC imaging studies are insuf¬ 
ficient for obtaining a full understanding of LyC emission 
from star-forming galaxies because they are so heavily 
contaminated by foreground objects. In order to elim¬ 
inate cases of foreground contamination, it is essential 
to obtain high-resolution observations of putative LyC- 
emitters to confirm the redshifts (either photometrically 
or spectroscopically) of all substructure associated with 
the galaxy. To date, all such observations have shown 
that candidate LyC-emitters from ground-based studies 
with anomalously high ratios of ionizing to non-ionizing 
radiation, both within our sample and in the literature, 
have proven to be from foreground contaminants. 

Future progress in understanding the physical proper¬ 
ties of LyC-emitters and the role of star-forming galaxies 
in reionization is contingent upon two factors. First, ob¬ 
servations of sufficient depth to probe ionizing radiation 
in galaxies at the faint end of the luminosity function 
must be obtained efficiently for a large sample of galax¬ 
ies. Second, these observations must be obtained at high 
spatial resolution, and with redshift information for each 
galaxy component. With such observations, it will be 
possible to distinguish between emission from foreground 
contaminants and genuine high-redshift LyC emitters, 
learn more about LyC photon escape by studying the 
multiwavelength properties of LyC emitters, and place 
more stringent constraints on the contribution of star¬ 
forming galaxies to the ionizing background. 
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Figure 16. 373 X 278 postage stamp images of the photometric LAE candidate lae^OlO. From left to right, images displayed 

include LRIS NB3420 (LyC emission), HST U 336 (a combination of LyC and non-ionizing LfV), HST Veoe (non-ionizing UV 
continunm), HST J 125 (optical, bluewards of the Balmer break), and HST Hiso (optical, redwards of the Balmer break). In 
the high resolution [/aae and Veoe images, lae4070 appears to be composed of two components (a and b, separated by 0724; 
labeled on the Veoe image). The photometry and EAZY SED fits for these two components are plotted below the images. In 
the left-most panel, the redshift is allowed to float during SED fitting, while in the middle panel the redshift is fixed to 2 = 2.85, 
the redshift corresponding to that of the redshift spike in the HS1549 field. The right-most panel shows the redshift probability 
distribution. Colors and symbols are as in FigurejO] The main qnalitative difference between the SEDs of lae4070a and lae4070b 
is the magnitude in J 125 . In the case of lae4070b, the resulting SED has the typical shape of a contaminant. However, in the 
case of lae4 070a, the resulting SED cannot be well-fit by model spectra at any redshift. As discussed in Section (5.3.21 it is thus 
possible that this is a case of a LyC-emitting LAE {lae4 070a) projected very close to a foreground contaminant {lae4 070b). 


of Hawaiian ancestry on whose sacred mountain we are 
privileged to be guests. Without their generous hospital¬ 
ity, most of the observations presented herein would not 
have been possible. 

APPENDIX 

A. OB,JECTS WITH AMBIGUOUS SEDS 

Here we discuss the interpretation of four SEDs of 
photometric LAE candidates with NB3420 detections 
where the SED shape is ambiguous, and there are no 
spectroscopic redshifts available to confirm that the ob¬ 
jects are indeed at z ~ 2.85. In total, we present three 
photometric LAEs from the Appendix of M13 {lae4070, 
lae5200, and lae6510) and one photometric LAE that 
was originally in the spectroscopic LAE sample, but for 
which the HST data showed that the wrong redshift 
had been assigned {lae7180). Eor lae5200, lae6510, and 
laellSO, Keck/LRIS spectroscopy with a total exposure 
time of 5400 seconds was attempted on clear nights with 
075—076 seeing, but no redshifts were measured. Eor 
lae4070, Keck/LRIS spectroscopy with a total exposure 
time of 17400 seconds was attempted under suboptimal 
conditions (intermittent clouds with seeing of 077—170 
during clear spells), but again, no redshift could be mea¬ 
sured. The lack of spectroscopic redshifts for these LAE 
candidates makes it impossible to confidently claim a 
LyC detection for any of them. Eurthermore, the shapes 
of their SEDs are ambiguous. Thus, while we can’t 


unequivocally confirm any of these candidates as LyC- 
emitters, we also cannot rule them out as foreground 
contaminants. 

The most promising LAE photometric candidate is 
lae4070. LAEs were selected by their H—NB4670 col¬ 
ors, as described in M13, and lae4 070 has a H—NB4670 
color of 0.82, slightly below the median value of the LAE 
sample and 0.22 mag above the selection threshold of 
H—NB4670>0.6. While lae4070 is fairly compact, it is 
still composed of two clumps separated by 0724 {lae4070a 
and lae4070h, indicated in Figure ITSl) for which we ana¬ 
lyzed separate SEDs. Both of these clumps are associated 
with NB3420 emission, and both clumps are detected in¬ 
dividually in (7336. Lae4 070b has the typical SED shape 
of a foreground contaminant. Lae4 070a, however, has 
an SED shape that could not be well described by any 
of the stellar population models we used with EAZY. 
The feature that most strongly indicates a redshift of 
z ^ 2.85 for lae4070a is the large break between J 125 and 
iJieo, combined with a flat Veoe ~ Ji 2 b color. In Figure 
[51 lae4070a has colors that place it in the same region of 
color-color space as typical LAEs and LBGs at z ^ 2.85 
(•7 i25 ~ Hiqo = 0.55, Veoe ~ >7i25 = 0.18). The emis¬ 
sion in the Usse filter is anomalously high compared to 
any model that provides a good fit Veoej "^ 125 ) and iLieo- 
While some of this emission may be contamination from 
nearby lae4070h, there is definitely emission in the U 33 Q 
filter at the location of lae4070a, which was not the case 
for any of the LAEs or LBGs without NB3420 detec- 
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Figure 17. HST Veoe imaging of lae7180, along with the 
corresponding Keck/LRIS spectrum. Thick red lines in the 
HST 14o 6 image indicate the location of the l”2 slit, and 
the LRIS 2D spectrum is aligned to match in the orientation 
and spatial scale of the imaging. As explained in Section 
15.3.21 the 2 = 2.93 galaxy near to lae7180 was not visible 
in the lower resolution LRIS imaging, and thus the bright 
Lya emission line in the spectrum was originally assigned 
to NB7180. However, analysis of the spatial distribution of 
the emission along the LRIS slit and how it corresponds to 
the HST imaging, along with the resulting SED fits of the 
sub-arcsecond components near lae7180, indicates that the 
neighboring galaxy is associated with the emission line and 
no emission line is visible at the location of lae7180 (which is 
now only a photometric LAE candidate). 


tions. If lae4070a is truly at z 2.85, then it is a LyC 
emitter. However, in this case, measuring the ratio of 
ionizing to non-ionizing flux for lae4070a in our NB3420 
image is impossible. The foreground galaxy lae4070b is 
located so close along the line of sight to lae4070a that 
there is no way to distinguish the NB3420 fluxes of these 
two objects with the 0'!7 seeing in the LRIS NB3420 im¬ 
age. Finally, we note that the foreground contaminant 
lae2436a (which spectroscopy proved to be at z = 2.04) 
has an SED similar in shape to that of lae4070a in Veoe, 
>1^125, and i?ieo, and this object lies near lae4070a in Fig¬ 
ure |^(Ti 25 — Hi6o — 0.60, Veoe ~ ,^125 = 0.06). The fact 
that degeneracies still exist in this area of J 125 — Hiqq 
and Veoe ~ >^125 color-color space where the majority of 
z ^ 2.85 LAEs and LBGs lie demonstrates the need for 
spectroscopic redshifts to resolve cases with ambiguous 
SEDs. 

We now consider the object lae7180, an LAE originally 
assigned a spectroscopic redshift of z = 2.930 based on 
a Keck/LRIS spectrum with a single high signal-to-noise 
emission line. The narrowband NB4670 filter used for 
LAE selection only probes redshifts of 2.80 < z < 2.88, 
so the redshift of z = 2.930 was anomalously high for 
a NB4670-selected LAE candidate. However, we con¬ 


sidered the possibility that a higher-redshift object with 
very large Lya equivalent width scattered into the LAE 
sample, and therefore retained lae7180 for analysis in 
M13. The HST imaging for this object, however, indi¬ 
cates that there is another faint galaxy in the vicinity 
of the z = 2.930 emission line, in addition to lae7180. 
Closer examination of the LRIS spectrum with respect 
to the HST images (see Figure ITTI) suggests that lae7180 
is offset by 1'.'13 from the location of the z = 2.930 emis¬ 
sion. With the HST imaging we were able to identify 
the true galaxy associated with the emission line, an ob¬ 
ject so faint (togos = 28.85) that it was undetected in the 
original LRIS imaging. The SED fit for this new object 
matches the redshift z = 2.930 identified in the spec¬ 
trum, and we conclude that the emission line belongs 
to this faint, nearby object and not to lae7180. Now 
that the z = 2.930 redshift is no longer associated with 
lae7180, this object can be reevaluated as a photometric 
LAE candidate with a NB3420 detection. 

Lae7180, lae5200, and lae6510, all exhibit ambiguous 
SED shapes. Rather than having SED shapes like that 
of lae4070, these objects have SED shapes similar to the 
one described in Section f5.ll which may represent either 
young, dusty, high-redshift galaxies or foreground con¬ 
taminants. As shown in Figure 0 for a spectroscopically 
confirmed galaxy with an ambiguous SED shape, EAZY 
gives a wide range of possible values for the redshift us¬ 
ing both the PEGASE and SMC-reddened BPASS mod¬ 
els. However, as Figure [5] also shows, the resulting red¬ 
shift probability distributions may differ when using dif¬ 
ferent sets of models. Thus, without spectroscopic red¬ 
shifts, it remains unclear for these LAE photometric can¬ 
didates with ambiguous SED shapes whether they are 
low-redshift foreground galaxies with an old stellar pop¬ 
ulation, or high-redshift LAEs with young, dusty stellar 
populations and LyC detections. While we cannot abso¬ 
lutely confirm if they are LyC emitters without spectro¬ 
scopic redshifts, the high-resolution HST imaging and 
SED fits can help narrow down the possible interpreta¬ 
tions for these objects. 

Figure [T51 shows the HST imaging and photometry for 
the sub-arcsecond component of each LAE that is asso¬ 
ciated with the NB3420 emission, along with the EAZY 
redshift probability distributions for these components. 
For each component, EAZY predicts a wide range of pos¬ 
sible redshifts (0 < Zphot < 4.5). Figure [TH] shows BPASS 
SED fits to the photometry both at low and high red¬ 
shift, and the fact that both redshifts can fit the data 
well demonstrates again the difficulties in distinguishing 
between low and high redshifts for galaxies with these 
SED shapes. In all of these cases, the detection of LyC 
emission depends on whether or not the LAE photomet¬ 
ric candidate is truly at z 2.85 - something that we 
cannot confirm for galaxies with ambiguous SED shapes 
and without spectroscopic redshifts. 

Here we consider possible causes for the scenario 
in which all three of these photometric LAE candi¬ 
dates {lae5200, lae6510, and lae7180) were incorrectly 
identified as LAEs. Lae6510 and lae7180 have fairly 
marginal E—NB4670 colors (0.68 and 0.63, respec¬ 
tively) when compared to the LAE selection threshold 
of E—NB4670>0.6. These values are in the lowest 15% 
of the E—NB4670 colors for spectroscopically confirmed 
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Figure 18. HST UizQV%o%Ji 25 color-composite imaging (5” x 5") and SED fits for the three photometric LAE candidates with 
the ambiguous SED shape described in Section [5.1 1 In all cases, the sub-arcsecond component shown is the main component 
associated with the LAE, and its position is indicated in the imaging. Results from EAZY are plotted in the three right hand 
panels. In the left-most panel, the redshift is allowed to float during SED fitting, while in the middle panel the redshift is fixed to 
z = 2.85, the redshift of the spike in the HS1549 field. The right-most panel shows the redshift probability distribution. Colors 
and symbols are as in Figure [G] Fits to these objects using the PEGASE templates yield slightly different redshift probability 
distributions, but they are qualitatively similar in that they span a wide redshift range (0 < z < 4.5) for all objects. 


LAEs in Ml3, and indicate that these objects are among 
the weaker LAE candidates and may have scattered into 
the LAE sample through photometric errors. As for 
lae5200, its E—NB4670 color is 3.21, an anomalously 
high value that may be due to contaminating light from 
a nearby bright star. Einally, these photometric LAE 
candidates were chosen for follow-up because of their 
NB3420 emission. As our new dataset reveals, true LyC 
detections in z ^ 3 galaxies are rare. Therefore, while 
many of the LAE photometric candidates from M13 may 
be true LAEs, singling out objects from the photometric 
sample with NB3420 detections may result in a higher- 
than-average selection of foreground contaminants. 

In summary, these four photometric LAE candidates 
with NB3420 detections all have ambiguous SED shapes 
that make it difficult to verify their redshifts photomet¬ 
rically and confirm their possible LyC detections. One 
object, lae^OlO, has similar J 125 — i7i60 and Veoe — J 125 
colors to many z ~ 2.85 galaxies in our sample and is the 
most promising photometric LAE candidate for true LyC 
emission. The other three objects {lae5200, laeSSlO, and 
lae7180 ) dis play the ambiguous SED shape described in 
Section 01 which may describe galaxies at many red- 
shifts. As we cannot unambiguously determine whether 
or not the four photometric LAE candidates discussed in 
this section are truly at z ~ 2.85, we adopt a conserva¬ 
tive approach and do not count the NB3420 detections 
for these objects as secure signatures of leaking LyC ra¬ 
diation. 

B. OBJECTS WITHOUT FULL HST COVERAGE 

Here we present postage stamp images of galaxies with 
NB3420 detections, but for which imaging was not avail¬ 
able in all four HST filters (see Eigure lTS)) . While insuf¬ 
ficient photometric data prevents us from fitting SEDs 


and determining photometric redshifts, we attempted to 
examine the morphologies of these objects in the Egos 
image, when available, to find objects with simple mor¬ 
phology where the possibility of contamination is low. 
However, none of the objects shown in Figure [19] have 
simple, compact morphologies in Ejoe- All objects either 
break into individual sub-arcsecond components or show 
extended diffuse emission. As both clumpy z ^ 3 galax¬ 
ies and foreground contaminants may be responsible for 
these multi-component Eeoe morphologies, we are unable 
to draw any conclusions about the amount of contami¬ 
nation in this galaxy sample. 
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Figure 19. 5” X 5” postage stamp images of galaxies with NB3420 detections, but for which imaging was not available in all 

four HST filters. From left to right, objects are displayed in LRIS NB3420 (LyC emission), LRIS NB4670—E (indicating Lya 
emission), LRIS V (non-ionizing UV continuum), and (when available) HST Uzm (a combination of ~ 80% LyC and ~ 20% 
non-ionizing UV) and HST Veoe (non-ionizing UV continuum). For objects with imaging in both HST Uzse and Veoe, we show a 
color-composite image of these two bands. Red (blue) circles (l”0 diameter) indicate the centroid of the NB3420 emission (Lya 
emission). The redshift of each object is indicated below the object name, or, if the object is a photometric LAE candidate, the 
U—NB4670 color is indicated. Postage stamps follow the conventional orientation, with north up and east to the left. As there 
is insufficient photometric data to ht photometric redshifts for these objects and as none of these objects have simple, compact 
morphology that would lessen the chance of foreground contamination, we cannot draw conclusions about the contamination 
rate of this sample of objects. 
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